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rORTWORD 
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TRW Space and Technology Group 

R. Avant 
Dr. F. Lee 
Dr. J. Tront 

Virginia Polytechnic Institute and State University 


Task 3 - Comprehensive Power System Analytical Modeling Approach 
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i:Xi:cinivi SUMMARY 

Powor suhnysToni analyT.lcal inodols havo ln'on list'd hy iinvci'niiii'iil and imlir.h’y 
for spacoiTafl. dositin l.o porforin i:onfi<iurarion lrad(’"OHV,, dot cnn inc pi'rlorinanci', 
spocKj " ’"iiial and oltH;lrU.aI inl t’rfai’.iv;, vt'i’ ify I i 1 y and del t>riiriin' clct - 

1 roimuinotic ooiiinatlbil ity niart|ins. fri'qui'itt. ly, inodols providt' analyl.ical data 
in support of 1n*orhit. anomaly rosolul.ion or prt'dici fuluri' capahilit.y consitii'r" 
iii(| dotp'adt'tl oporal.ion. 

Modi'Is will play an important. rol(> in ilosi»)n and vori < icat.ion of ful.urt' 
lar^o space platform and space station power subsystems. It is hert* that the 
need is most critical because design verification throutih integrated systems 
tests will be prohibitively expensive or impossible. 

Consider a large space station with an orbital lifetime of 10 to 30 years. 
Experiment and power source modules having power characteristics unknown at 
the time of space station power system design will be launched into orbit by 
the srs and attached to the orbiting station. It is planned that new modules 
be attached to the station many years after station launch. Discovering power 
system incompatibility at the time of attachment will be expensivi* since an 
STS flight must be repeated if in-orbit repair is not possible. Instability 
of the power bus due to incompatibility may result in a dangerous condition 
for the entire station. 

Therefore, a comprehensive power system model which is more accurate and 
flexible than today's models is required to verify through analysis, rather 
than integrated system tests, large power system performance and stability. 

This is a significant departure from today's integrated system verification 
through testing prior to launcl.. This report addresses improvements needed 
in current models and provides a roadmap to development of a compiM'hensive 
power system model. 

SCOI'i: 

this study examines spaiecraft power sub';ystem models as doeument('d in the 
literature and as determined from a survey of government agencies and industry. 

I he st ate-ot ' the-art is compari'd to desiri'd compri'hensi ve powi'r siil> .ystem modid 
attributes. Weaknessi's or inadt'guacies in current models are identified as 
art'a.. tor improvement. An iippro.uh to comprehensi ve power subsystem model dc'velitp 
ment is pri'sented togt'ther with n'commeiul.i t ions lor test d.ita b.ise dt'velopment 
re(|uired to verity and validate the model. 
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OBJECTIVES 

The study ohjoctivos wore to dotormlno the capabllltlo;; of power system 
modeling techniques presently used by government and industry; to characteri/^e 
available AC modeling capabilities for power system comporionis typically 
utili;;ed in 2-15 kW power systems; to delineate component modeling technique 
improvements, permitting accurate simulation of AC and DC characteristics 
of power system components; to develop an approach for orderly ac( onipl ishment 
of necessary analyses and device testing which support developmeni of a 
comprehensive power system analytical tool. 

SUMMARY OF RESULTS 

We recommend a set of four fundamental model types, each of which 
performs a different, essential task. Together this program set comprises 
a comprehensive power subsystem model. 

1. A power subsystem sizing and synthesis program, capable of 
estimating cost, mass, volume, area, and other attributes 
of a single-point design. This would be used only during 
the conceptual design phase of a spacecraft program. 

2. A DC model of the power subsystem and its interfacing sub- 
systems. This is used during phases* B, C, D, and E'of the 
spacecraft program for studies of power consumption, responses 
of the subsystem to environmental changes, and prediction of 
steady-state voltage and current throughout the subsystem. 

3. A small-signal AC model of the power subsystem and its com- 
ponf *■‘ 1 . This is used in phases C, D, and E for the purposes 
of esi..matim subsystem and intercomponent stability, bus 
impedance, and for determination of responses to small-signal 
transients . 

4. A large-signal transient model used during phases C, D, and E 
of the program for the purpose of determining the response of 
the subsystem to large transients such as state changes and 
faul ts . 


1 
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In the preceding paragraphs, model requirements have been defined. In the 
following are defined certain essiMitial attributes or characteristics which 
all of the models must have in common to create a comprehensive modeling 
and analysis capability. 


★ 

The mission phases are defined in Section 3.0 
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a) Coinninnal it y ami cuiii(Mt iliil i ( y: I ,i< h of flic lyjji's of moilols 

iiuir.f iiuilnfain a i.oiiiiiional ify of roforomc with all of fho ofhors ; 
i.o», tlioy lUMsi fill prudli I. l:h(> sniiu' iiH rlhiitos (ir)d porroriiiiifu p 
for flu? r.aiiio ‘iulu'.y'ifoiii. 

h) Motlularlly: I a( h ol I ho i.oiii(iOiioiif:> of fho pow(?r r,iih!iyi;f(!iii 

iihoiilil ho imulolod as <in 1iulo|uuulotif iiiodiilo havimi tho i’ol lowlmi 
affrihuto!;: 

• Capahl(? of h<?1iu,j o|u?rafod iiidopoiuloiil.ly a;; a coiiipoiuuil iiiodol . 

• Capahio of hoiiig infc'^rafod Info a powor nuh'.yiifoin iiiodol . 

• Capahio of rcplacltui or lu?ing roplacod hy an altornalivo 
model havinq differenf input requirements. 

• Data ha;:e - each module must be provided with an independent 
data base specifically suited to its own needs. 

c) Efficient use of both core memory and computation time. 

d) Verifiability: Models should be verifiable as independent 

modules and as a complete integrated power subsystem. Veri- 
fication data base must be provided independent of the model 
input data base. 

e) Operational simplicity: The models must be designed for use 

by a working power subsystem engineer whose understanding of 
tho program and computational facility are limited. Clearly 
stated, complete user documentation is essential to the 
successful use of such a complex set of models. 

Consideration will also be given to the development ot a program control 
format based upon MENU selection of program functions. This is desirable on 
systems containing CRT displays, and the potential advantages of such a 
system in terms ot user simplicity and compatibility are great enough to 
deserve consideration. 

Models may be divided into two categories: existing models which require 

only adaptation to the program structure, and new or improved models which 
require development. New model development proceeds through trade studies 
to select the approach to be used, development of the mathematical or logical 
algorithm, and coding. Irade studies and algorithm development can bi> done 
prior to or concurrently with the program structure development. Coding 
must aw.nt the dt*velo|)ment of program and data bas(> structure before it can 
be accoiiipl ish(>d efficiently. 
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Tor thoso models whovio concoRt'.; aro woll iindnrotood, a data ai iiuisIMon 
pt’ocirain can bo Initiatod indopond('iiM y of l.ho d(*volnf)iii(*n1 of tho final iiiodol 
codo. Howovor, whom tho ntodol alnorlttiiiir, havo yoi io bo dovolopod, iho form 
of tho data roqnirod, and in noino casoi;, tho kind of data ro<|uirod to providf’ 
an InfMft to Iho modol aro unknown or 111-dofinod. for thoso ca<ios. It will 
bo nocosi’.ary to dofor tho dnvo1o|»inont of a data acquisition proqrani until 
tho alqorlthni Io bo usod has boon dofinod. 

fiquro 3-3 (paqo 3-13) shows a suqgostod schodulo for modol atid data baso 
dt.'volopmont. This schedule is flt?x1blo, and Is capable of respondintj to varia- 
tions In funding availability, NASA priority of interest and other factors. 
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1.0 INTRODUCTION 

Tho qrowtti In tho and coiiiploxITy of spaccrr’.iTI- power •■vsteiiis , 
coupled with hl<|lior o(|uipiiient '.wUcIvIno Trecpiem In-, and ati 1ncr(M'.e In 
payload sensitivity to hus noh.e, has focused attention on a current major 
deficiency - the ability to design and test latsp! power systems. Ix1st1rif| 
analytical models of spacecraft pctwer systems have not Kept abreast with 
these evolviny rei)u1remnnts and, as systcmis orow In sixe, there are 
practical limits In the ability to (jround t(!St fully Intefjrated power 
systems. Roccfit anoiiidlous behavior experienced on both NASA and Air rorcc 
satellites underscores the need for better model Itu) techniques lor both 
the design of stable power systems and the efficient management of power 
during tho mission. The development of accurate analytical modeling 
techniques was cited as a highest priority Item by the Power Subsystems 
Panel at the OSTA/OAST Flight Technology Improvement Workshop in 1979. 

Uncertainties in the analytical modeling of power systems are 
derived from two basic sources: first, the component characteristics, 

especially the AC characteristics, are unknown or at best poorly defined. 
Consequently, DC models are used for most analyses, coupled with approxi- 
mations for solar array simulations. Secondly, all-up systems testing 
with illuminated solar arrays has become impractical for mul t1 kilowatt 
systems due to the large area of Illumination required and the risk of 
damage to the fragile lightweight structure deployed in a gravitational 
field. These problems will be significantly aggravated for the even 
larger systems projected for the future. 

In those cases in which the inadequacy of DC power system modeling 
has been recognized, dynamic analyses have been performed. However, these 
analyses have proven to be both time consuming and e::pensive. In addition, 
the analyses have been handicapped by unknown or uncertain component 
characteristics. Only through a more thorough knowledge of both the AC 
and DC characteristics of the devices which make up tho power system com- 
ponents can the AC/DC characteristics of the components bo accurately 
simulated. With the accurate modeling of the components, future power 
system design and evaluation can he accomplished by the synthesis of the 
analytical system model. The development of the analytical model for 
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powor syatoinr. will ulr.o providr^ for accuratn r.liiiulrition of tho solar array 
and battory, porin1tt1ri<i roallsflo nround fosliiKi of tho systoin and rolifiblo 
sytifoin dovolopiiiont. This dovolopiiiont would f.icllitato l.ho analysis of 
sovM'ral approat.hos which could satisiy a sot of uivoii rooulromoni s, and 
allow optimisation of tlio ovt'rall powor sy'.loiii. 

Tlio study was dividod into thr('o basic I, asks which aro fku union i.od In 
this final report : 

'* ) T‘'jsK. J Power, System Model liiy Teclinl < tuo'; 

Conduct a stato-of-tho-art survey combined with a comprehensive 
review of power system computer modeling techniques/approaches used 
by Industry for performing an accurate simulation of system perfor- 
mance of both earth-orbiting and Interplanetary probe spacecraft. 
Identify the significant capabilities and drawbacks of each analytical 
technique, along with the areas in modeling which require significant 
improvements. Determine the simulation adequacy of dynamic load 
changes on the power bus performance. 

Task 2 - Power System Component Modeling 

Determine the adequacy of each AC and DC characteristic modeling 
technique for the components typically used in power systems operating 
in the 2 to 15 kW range. The improvements required in each of the 
modeling techniques will be specified, along with the degree of un- 
certainty associated with the present modeling approach. 

c ) Task 3 - Comprehensive Power System Analytical Modeling Approach 

Outline the necessary procedures to permit the development of a compre- 
hensive power system analytical model. This will occur after the review 
of the available analytical modeling techniques and will also include 
the areas in modeling which will require upgrading. Also to be identi- 
fied are any power system component or device testing required to 
assemble the AC and DC characterization data base from which the 
component and power system models can be derived. The recommended 
or suggested order in which testing and analysis should bo performed 
will be provided. 
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1.1 St.iuly OhJocTiyo'', (HkI 

ITin prliiiiiry nb.|ot:ilvo of this 1iill:1«l of tlio powor '.ultsystniii 

inofinlinn f-lnvnlofiiiioni: lr, t.o cloflito tho roqLriroinpntn of (-oiii|irf»lioii'.lvo 
iiioclol (or f.ni, of iiiodol'.) for solar array haffory olnctrU fiowor -.ultsyfitoin. 
To ac.c,nnip11'.h fills fho slu'ly was dlvldod liif.o tlio foriowliKi tasks: 

Task 1, DotoriiiliiP tlio ( apabll It loii (and dcfltloiiolos) of tlio I'owor 
suhsyfitoiiis iiuulol liKj tochn'liinos iirosontly 'ii ir.o. 

Task 2. Dotoriii'lno tlio liii|»rovoiiiunt ro(|irlrod of (*aoli iiiudol . 

Task 3. Doflrio tlio ronulrcnicnts of a coinploto power subsys tom 
iiiodol , and tho procoduros for dovolopmont of a iiiodol 
and data baso. 

1 2 S tudy Resul ts 

1.2.1 T ask si Jt 2 “ Pow er Subsy stem Models Present 1 7 In Use 

Recfii 1 red Impr o vements 

1 . 2 . 1 . 1 Steady State DC Performance Models 

It is in this class of complete power subsystem models that the most 
extensive body of work was encountered. Most models are similar in nature, 
consisting of the major power subsystem components modeled empirically and 
tied together by a spacecraft-unique model of power subsystem logic. They 
have no AC or electrical transient capability, although some are capable of 
simulating thermal transients. 

1 . 2 . 1 . 2 Small Signal (Linearized) AC Models 

Such models are used primarily for stability analysis. Individual 
component Thevenin-equivalent linearized models are available for most 
components, but vary considerably in validity and applicability. Combi- 
nations of two or more such models have been made for specific stability 
analyses. All-up power subsystem models of this type have not been found. 

1 . 2 . 1 . 3 Tran sient iNon-Li near) Model s 

Transient models of adequate accuracy are available for only a few 
power subsystem components. An all-up transient model of a direct energy 
transfer power subsystem has been simulated on a digital computer (Reference 
67). 
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1 

1 

I’dWt'i' (dill Co(ii|ioiii'iil. MiKldl'i ; 

R.il I nr i(". . 'idvoi'dl iiimlcls iit hdl l.cry (H’rliinii.iiiM’ Iwvi' i>n( oiini rrcfl 
1ii tlid M li'i’.il iiro. 

1 . I (ml'.' ii|i I .iltlc IlHidcl (A' I'M)* I , 

I ' 

I*. I (|ii ivd li'iil I In. Ill I iiimli'l (A' lV) I 

I lllp lr l( d l WdVO'.lldpn IIKIlll'l 

4. IIOMSAI/Ill I Inrlidck inudtil (A-I'l’) i 

Hid 1iu)k»up Idbld iiimidl Is d DC st.cddy“n l.(il,d iiioildl ( t.lii’riiiii I I i dir, Idtii.s cdii I 

l)d fo'Muwdd), diid lias hucn widuly used by many companlys. The iiiodol tidods |, 

liiiprovdiiidiit. in dt'i icidiicy inoddlimj and lidai (jdiidralion rate comt. '.ai lo!. 1 

I 

The models are not universany applicable. Diffen'iil data ba<- ■ r.'jst j 

be used for different orbits, and to simulate battery wp'i t. ii.e ! 

models are available only for NiCd batteries. Ma s .iiodel of NiH., 

cl 

or Apll., was found. 

u 

The equivalent-circuit model is also a steady-state model. It has 
the advantage of greatly reduced data storage requirements. As far as can 
be determined, this model has not been pursued since being reported in 1970. 

The empirical waveshape model also has not been pursued after initial 
development. The pressure predictions of the model are not accurate and 
require improvement. It appears to be valid only for low earth orbit (LEO) 
appl ications . 

The COMSAT/Bil Icrbeck model (Reference 22 ) considers only discharge 
and therefore would liave to be expanded if it is to be included in an 
overall power system model. 

Rower Distribution 

Two modeling approaches were encountered: 

1. Opl imi.;:a tion based on mass 
Optimisation based on cost 




*A-XX ri'fers 


to Appendix A document 


referem ('. 
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In the first caso, the existing models do not account for fuse, switch and 

connector losses, which In some cases could override conductor losses. In 

the second case, the optimum conductor cress sectional area Is determined 

as 0 function of the cost of the conductor, energy storage, solar array, i 

thermal control and power conversion. Both models can bo adapted for use \ 

by an overall systems model. Neither 1s suitable for use as a performance ; 

model, el tiler Ad or DC. 

Solar Arr ay 

Because of the Increasing sUe and reduced specific mass of future : 

solar arrays, all-up testing may not be feasible, and accurate modeling 
becomes Increasingly Important. i 

Existing steady-state DC models of solar array electrical performance 
are generally very good, although shadow effect treatment and multiple sun i 

Illumination levels both require Improvement. ) 

■ i 

Small-signal AC models are available, giving good first-order results. ! | 

They are relatively easy to use. The needed data can be measured or computed | 

without great difficulty. Considerable improvement in model results should ' I 

i 

be achievable by the addition of parasitic wiring capacitance, capacitance, i 

cell-to-substrate capacitance, and photo-inductance effects to the model. I 

I 

Three general approaches were found; | 

■ i 

1. Model each cell Individually. H 

\ ; 

2. Model the array as a multiple of a single cell. 

3. Combinations of 1 and 2. 

The first approach found provides a highly accurate performance analysis 

of an existing solar array, but requires voluminous data input and exces- | 

sive measurement and computation cost. ’i 

I 

The second presumes a single, uniform cell characteristic. Because !| 

of the largo fiumbor of cells and their well controlled characteristics, |.j 

this approach models the characteristics of the array (without failures) |{ 

nearly as well as the first. However, its handling of an array (ontaining 1| 

failed or shadowed cells is inferior to that of the first. j; 


No large-signal transient models were found. 
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iu;-(l(' rmivot l I'r 


lltroo appro.u lt(". to ( diivt'i’l i>r iikuIpI imi witp ith'iil i I ((••I ; 
I . ni'.i ri'l I' ttvt'i’.nh' 


. Hi stTi'i(' ( inic ilom,i i m 
. 1. CuMItMMl t'M‘t:ilit .1 Mill vs is 


Ihc ,llsi:ivr,. h. i„ „M', ,iml „p hair 

or tho swir.hio,, fro.|„o.u,v toi- /U: analvvl.,, 1 Imo ,lomai„ mo, I, 'I 

K a .l,■la 1 l,■.M„r,,.,■nal woi-kl,,,, , „s,,r„l n„al ,lo..l,|o ,:h,., k. II,,. 

.|,',i<.,.al a„aly,,ls a|,|„„ad, is v,-,-s.nri,> .„„l ,.asil,v oxpaodal.l,., how- 

ovin* tli(' cost iiuTo.tsos wil li ifUT(\isino foinpl I'xi | y. 


KtHiul.dors 

StiMdy-stdt.e DC niodol s havo tn^on for two typos of shunt s, tho. 
|Mrl.i.n-lino..r, .oul tho full soMuontial stuiinontod linoar shunt, usod in 
IjSCS 1 1 , tuul III AO, rospoi 1 iv('l y. 

An analytiial iiuulol lor AC, small signal analysis docs exist for tho 
shunt rooulator. Ihis "shunt loop c.ain modol" has tho advantages of simpli- 
city in final term, usotulnoss in predict iiu] shunt suhsystem stability, and 
ease of application to computer programs which accept s-domain transfer 
tunctions or differential equation inputs. It has, however, substantial 
limitations, some of which are as follows; 

* model required to establish the 


• It is ilitlicult to alter I lu' model. 

• !apab'iHly* Provide larqo siqnal transient prediction 

• No powi'r dissipation or elt'im'iit lailurt' analysis is provikled 

• It is .ipplicabit' only in I lu' limsir t>peralional imuh-. 

• No uvntrol output pissliiMrui is pisivirlod. 
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To provirto ..do.|Ual.o «iodol1n,| 1 1 ty for fho simnt (or parttol sliuot) 

roouktor lo a -.alollUo powrr poorroflon ■oih..y-.tom osinp a Solar /Irray 
SwUrhiloi UoH (SUSII), a network approarli to modnl1nt| may lio nioi r appio 
priato. fhr «|olvalrnl nol.work ap|o-oarh lor. I ho advaiitagos of rr.pilrlmi 
al ioUtal offort. is oastly aUorod, providos easy 
individual .ompononl behavior, and allows laroo sltinal. r.mall s1.|oal ,ind 
DC porformanoo analysis, as well as power dlsslpallon predictions. Tho 
network approach has dlsadvantapos , In that a sophisticated computor propram 
with large network capacity Is required, and that parameters tor Individual 

element models must be obtained. 

Building blocks exist for tioth approaches to shunt regulator modeling 
(1 e element models, element s-domain descriptions, and analysis programs. 
The synthesis of these blocks Into an effective shunt regulator model to pro- 
vide adequate computer analysis capability has not been accompl.shed at 

this time. 


Power Control Models 

The Individual component models described above can be Integrated 
Into a complete power subsystem through the use of a power-control model. 
The spacecraft-specific power control model embodies all of the power su 
system control loyic, including: 

• Battery charge/discharge controls 


• Switching controls 

• Eclipse - daylight controls 

• Sonsor elomonts necessary to control subsystem logic. 

Most of the power control models encountered have had only DC capa 

bilily. 
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1.2.2 Task 3 - MqOol Dovq.1 ohul'H t. Prococliinis 

For the develoiwenf of an neenrnto power suhnystc* model, the following 
are roquirod; 

• A inodol of each power subsystem component, coniplete with model 
algorithm, and an adequate data base. 

• A floxible dynamic program structure capable of accommodating 
5 Tide variotror component models, some of which are presently 
ill-defined. 

The eemprehensive model is envisioned as a set of compatible specific: 
purpose programs, each of which calls upon a library of power subsystem 
component models and data to model a wide range of power subsystem configu.a- 
tions. This approach minimizes both complexity and cost, while covering e 
entire range of model requirements, from the simplest sizing model through 
13C and AC performance, to the most complex transient model. 

Driver program structure will require a flexible means of data storage 
and communication, so as to permit use of a wide variety of component models. 
Initial effort will focus on program structure development. This will 
permit the program to be used as soon as the first component models are 
available. This ensures program availability in the third year of the 

five-year effort. 

Test Contpq n^n t s 

Component testing Is required to generate an adequate data base. 
Battery testing Is required to Improve the efflclency/state of charge mode 
and the heat generatfon/pressure model; and to expand Its usefulness to a 
wider range of orbits, voltages, temperatures, etc. 

In the area of power electronics, testing Is required to clmracterlze 
devices based on «, urging technology, such as high voltage hybrid switches 
to be used for solar array reconfiguration and power bus control, MOSFETS 
to be used for digital .:ontro1 of solar arrays, and Integrated circuits to 

be uiicd to control switches. 
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Solar array testlnq is rnquirocl to develop an accurate transient 
model. The dynamic impedance is a function of voltane and currc'nt. A 
bettor underutandinq of plasma effects is needed before a model can be 
developed. 

Analysis or testing is required to improve the discrete average 
model of the duty cycle modulator. The inaccuracy is particularly pro- 
nounced for high performance converters employing multiloop control. 

The model is also unable to handle abnormal converter operations such as 
peak current protection mode of operation, saturation of magnetic compo- 
nents, and saturation of op-amp in the feedback control loops. 

The discrete time domain model is relatively ineffective for large 
scale system simulation. The general circuit analysis technique is 
fairly well developed in SPICE2 and SCEPTRE which are available computer 
programs. 
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2.0 POWTR SYSTKM MOOTl 1N6 TCCIINIOUfS (TASK 1) 

Fifluro 2-1 sliows tho work flow for Idsk 1. A UForal.uro soardi wa;i 
p(?rforiiif!(1 In ordor to rovlow oxistini) powor systoin coinputor iiiodolinq toch- 
niquos usod by industry to simulate spacecraft powor system porformance. 

Over 200 pnpors/documonts wore Identified arid cateqorized into one of 
five areas: 

1. Powor Supply Electronics 

2. Batteries 

3. Solar Arrays 

4. Solar Array Switching 

5. Power Systems 

Each document/paper was reviewed, and those having power system modeling 
applications listed in Appendix A. For each reference, the significant 
capabilities and drawbacks of the analytical technique were identified 
together with areas which require significant improvement. This informa- 
tion is summarized on technology evaluation summary charts. Each chart 
includes source identification, purpose and model description, capabilities, 
limitations and constraints, and areas of improvement. Technology 
evaluation charts are in Appendix A. 

In addition to the literature search, an industry survey was conducted 
to determine on-going effort in power system modeling techniques. Results 
of this survey are in Appendix C. The following sections present a summary 
of the findings. 

^ ^ ing (Figure 2-2 ) 

2.1.1 Battery MqdeJ s 

Several battery performance models have boon encountered in the 
1 i teraturo : 

1. A 'ook-up table model proposed by Bauer (A-28) in wiiich the 
characterist its ot the battery are stored in a largo data bank. 
Liiu'ar intc’rpolat ion and doublo-inti'rpolation routines are used 
to find battery voltages, efficiet\cies, etc. The batt»'ry model 
is integrated into a complete power subsystem model by use of 
computerized pseudographical methods. 


r 
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An oquivalont circuit, inodol proposoil h.y 7 liniiiorninn atul t'Hrrr.on 
(A-27) in which 1 lu' hattory 1s rnprcsont o 1 h.y an oi|u1valnnl 
clrcull of capacitors, (lioilcs, resistors, atui switches. 

An tinpiihl Islu'tl iiiodi’l proposed liy Hauer in which I hi' current - 
vo1ta<)o characteristics of the liatierles are .v'prosented hy 
a set of empirical ('quatlons designed to reproduce the shape 
of the hatter.y charqe and discharge curves. 


4. A model of hatlery dlschanje vollatje .is a function of i;yc1e 
life and temperature Is reported hy hillerheck (A-l’2). 

Table 2-1 details present model capahilltlos with respect to several 
attributes, and Table ?.-? describes existing models uncertainty/ Inadequacy. 

Look-up Table Model 


This model consists of several independent models, each of which 
is designed to represent a different characteristic of the battery cell. 
These are: 

1. Data bank consisting of battery voltage as a function of 
relative current [current/cell capacity (A-Hr)], temperature, 
relative state of charge (percent remaining charge), and state 
of charge at the end of the preceding discharge (relative). 

By entering this data bank with temperature, state of charge, 
and depth of discharge arguments, and performing multiple 
interpolations, a single current-voltage curve for the cell 
at the desired set of conditions is obtained. 

2. An efficiency model consisting of the incremental cell 
efficiency as a function of relative charge current and temp- 
erature. This table is valid only for the case in which the 
battery is charged from zero state of charge, and must be used 
with an accompanying algorithm to determine efficiencies during 
recharge from higher states of charge. 

3. A battery heat generation model. This has taken several forms 
as the programs evolved, and is now an approximation oi the 
thermodynamic heat generation properties of the battery cell, 
including energy storage as oxygen pressure. This permits 
modeling of' the delay in heat evolution experienced in actual 
battery cells on repetitive cycling. 
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• JIC^E^S FOR 




• ^jCOELS FOR 


• ACCEPTABLE ACCURACY FOR GEO ORBIT 

« 

• Heat generation and temperaitjre mcdels exist 

-- TESTING/ANALYSIS/YERIFICATION required for multiorbit CAPABILITY 

— TESTING/ANALYSIS/YERIFICATION required to OBTAIN ACCEPTABLE 
PRESSURE MODEL 

• CELL AGING MODELS AVAILABLE, TESTING ReXlIRED TO IMPROVE ACCURACY 

— FUNCTION CYCLES^ TEMPERATURE 

• NO KNOWN FAILURE OR TRANSIENT MODELS 

— testing/analysis/v’erification required 

• MODELS AVAILABLE^ THEVENIN EQUIVALENT 


Note: above coftiENTs apply to NiCd baheries, all <\reas 

REQUIRE testing TO ACQUIRE AN ADEQUATE DATA BASE 
FOR 01T€R. TYPE BAHERIES 


Table 2-1. Component Model Capability Summary 

Batteries 
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TABLE 2-2 

Battery EIodel Uncertainty/Inadequacy 


• HEAT GENERATION - CELL HEAT GENERATION' MODELS WEAK - OFTEN GIVE POOR 

RESULTS WHEN ORBIT PARAMETERS CHANGED 

• MULTIPLE ORBITS - MODELS FREQUENTLY DO NOT GIVE STABLE RESULTS AND 

DRIFT, RESULTS NOT IN AGREEMENT WITH FLIGHT DATA 


• CEa TYPE - NO COMPLETE NjH^ OR AgH2 DATA BASE EXISTS 

• RECONDITIONED CELL DATA BASE 


• BATTERY CELL FAILURE MODEL/LOAD SHARING SIMULATION 

• pu .SE discharge/charge data base 


Tho 1ook"U^ ntodol har» boon loiod aftor varyinq donroos of modi f icatluri'i 
by snvoral Institutions, with a varylno dooroo of suocosf.. ,)PI has niodiflod 
It for usn with Intorplanetary r.pacocraft. An Indopondnntly qonoratnd 
tabular modol war. mado for tho SKYLAB powor subr.yr,tom subr,t1tut1n<i an avorafio 
rocharqo ratio officloncy niodol for th(* Incroniontal iiiodol of coll offlclonclor,. 
NASA/6SFC usos a similar look-^iip modol for qonoral powor subsystoiii ovaluatlons. 

Mpdpl Capabllltl Whon coiiiblnod with a battery charqe model, a 
thermal model, and a power source model, the battery model is capable 
of making mlnuto-by-mlnute predictions of battery voltage, current, 
temperature, state of charge and, potentially, Internal cell oxygen 
partial pressure. In effect. It "flies" the battery through Its 
mission, predicting all of the Important characteristics of the 
battery In relationship with Its environment. With additional 
programming. It Is capable of simulating the performance of a battery 
or batteries containing one or more failed (or partially failed) 
cells. 

Model Limitations ; The look-up table model has a numbev> of limita- 
tions: 1) It Is Incapable of modeling electrical transients. It 

Is useful only at steady state. (Thermal transients, except in 
masses of very small time constant, can be followed accurately, pro- 
vided that the interval between computations is short); 2) It does 
not account for actual capacity degradation due to loss of active 
materials, or conversion of active materials to the inactive state 
with cell usage. Instead, these losses are modeled as a voltage 
degradation, and require replacement of the current-voltage data 
set with a "degraded" set. The Input data for the characteristics 
of the battery have to be modified for each set of battery conditions. 
Since these data consist of voluminous five-dimensional tables, It 
becomes expensive to generate a lifetime or cycle variant data base. 
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Accuracy connidorati onr.; Exporlnncn at TRW hac Ind to tho following 
ontlmator, of accuracy. 


Qattory voltago 
Battory current 
State of ehargo 
Moat (jonoratlon rate 


l.fj'^ 

l.!V,: 

5.0 » 10 . 0 ';:; 

Variable, doponding upon 
•jtago of battory chanjo 
or dischargo 


Potential for Improvement : Areas in which liiiproveinents could onhance 
tho usofulnoss or rango of applicability of tho progranir. are: 


• Development of an algorithm for manipulation of the voltage and 
efficiency data bases to reflect the variation in cell charac- 
teristics with life or cycle histovy. This would Increase the 
range of applicability of the program. 

• Development of an Improved efficiency model. This would 
enhance accuracy of heat generation results, and state of 
charge results. 

• Development of an Improved heat generation model. 

Equivalent Circuit Model of 7.1iniii ernia n and Peterson (A -?7 ) 

This model consists of a battery steady-state equivalent circuit 
comprised of two or more large capacitors, a bidirectional diode pair, and 
a network of resistors. The bidirectional diode pair simulates battery 
cell voltage hysteresis and battery heat generation. 

jt.i_^.^ I The model simulates the general shape of the 
battery charge and discharge curves, and the offset between them. With the 
addition of considerable complexity, it appisirs That the same mod(*l might 
be designed to simulate both steady-state and small signal AC behavior and 
heat generation hysteresis caused hy nickel -cudmium tell oxygen storage. 

The amount of data storage reguired is small. 
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Thn orlflinrtl popor rtpp(^(4rr)d In 1470, nud nn furrhor dovolopnmnl. hds Oonn 
roportod Imidinq to thn tuncUi-ilon thrit 1n npltn of Itr. prniiilr.o, thn work 
Wrtf. not purr, Mod furthnr, or wnn unrucrnnr.ful . No coinpnr lion 1i ihown bn- 
twnon fl 'lflht or toil, dnto (uid thn i.ninputnr prod lot Ions , 

Moclol J Onn of thn d1 fficiil t ini proK-nl.nd hy thi'. inodol li 
tho roqulrnimmt for Input dntn 1n u forin which li not dirootly nioftfiuruM&. 
fho dotti oonvonlnn noniploxitloi iiro lorrin, niid ditto /lotiuli I Lion unwlnldy, 

Ar, fliodol liiiprovonionti am mado tor miiari AC nignal por forma noo. for oxygon 

Rtorago iiiodoling, oti;,, thoio coinp1cx1t1oi art- Mknly to Incmaio gooinotrl- 
cally. 

WtLyt-lh'iijj.n Model (A-?H) 

The "Empirical Waveshape Mode!" assumes the battery charge/dlscharge charac- 
teristic wave to consist of three regions, a sloped plateau line in the 
mid-region, and two exponential functions added to the plateau to simulate 
the beginning and end of charge or discharge. The combination of these 
three functions model the zero-current cell voltage as a function of its 
state of charge. This "zero-current" voltage is then compensated for charge 
or discharge currents by adding or subtracting a voltage Increment whose 
value is a function of the relative current. 

The overall model contains a battery pressure and heat generation algorithm 
based upon thermodynamic relationships, and a tabular efficiency model 
similar to that described in the look-up table model. 

Mo de 1 ca pa b n i t ji es ; The model will predict battery voltage, current, 
heat generation rate, state of charge, and internal pressure of the 
battery cells. It is capable of expansion to include the modeling 
of partially failed cells. 



H 

i 
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Model limitations ; The model requires an indirect and In some cases 
a trial -and-error derivation, of Input parameters from available cyclinq 
data. In testinn, the model maintained qood aqreement with actual 
cyclinq data over the ranqe of 0 to 26"C in a ninqle low earth oH.!i 
cycle comparison with test data. It has never been compared with tost 
or fllqht data at other temperatures, or at varyinq depths of discharge. 

In multiple cycling runs It may require as many as 30 orbits for the 
model to converge from the assumed set of conditions to the final con- 
ditions In which equilibrium Is achieved. Similarly, a large number of 
orbits may be required before It settles down as a result of a large 
change In system loads. This Is not unlike the behavior of actual 
spacecraft. 

teura cy cons 1 derat 1 on s ; In the single test In which comparisons were 
made with low earth orbit test cycles, voltage, current, state-of-charge 
agreed with test data within 1-5%. Pressure predictions In a few test 
cycles showed a pattern of pressure variation in phase with that In the 
test cells but different in magnitude by as much as 50 to 100%. 

Potential for improvement ; The model fails to take into account the 
variation in current density with state of charge at constant current. 

The assumptions regarding the rates of reaction associated with the 
pressure and recombination reactions are oversimplified for best accuracy; 
further improvement is possible. 

COMSAT/BI llerbeck Model (A-;’;’) 

In this case the battery is modeled as a voltage generator whose output 
Is a function of depth of discharge. In series with a fixed resistor, and 
compensated by temperature and aging factors. The model considers only battery 
discharge, and must be calibrated to the power subsystem prior to 
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MihIi'I i i',(»dhi I M it's : (liven nn .I'.'.itineil In. id t iiiTeiil, .in tis-.iiined ‘.|mi e- 

n’.s 1 1, d(’'.it|ii, .uul field oik'imI. inn d.iln Iroiii llie '.|MC('n’.ii I , tlie model 
will in’edirl li.itlery voll.n,)e .il end of dlMli.U'iie (iiiiniiiumi h.iMery 
vnll.ii|e . 1 ', .1 Inin I Ion ol Ihe nniiihet’ ol l yrli". in oi’Mt, .ind h.iMery 
e.ru' 1 nm|U'i'.i I to’e) , An dddition.il nliior I Him is ii'.od fo iinlude Hie 
('fftn.l,!, of liaftery nH’.ondU ioninq. 

Modi'l 1 iiiri fill ions : lleciiiise flu* model doi's nol eonsidei* battery i liarqe, 

11 will not pri'dii l pi'i'foniianee undee varyiini idi.irqe eonit'ol i ondi I inns 
nor will if follow Lliormal l.ransii'iil s. No heal qi'iieration data or 
jiressiiri* prodicfions are made. The iiunlel will nof follow AC variafions. 

Model aeeiirae.y: The reporfed model aeeur.uy ajipoars fo vary hefwoen 

;V;. early in lil'e fo 10-15’". later in life, wtii'n ealihrafed aqainst. 
early life data, and not includiiuj reeonditioniiiq. After several 
years of operational data have been acquired the model can he matched 
very closely to these data, and extended predictions made from this 
point . 

1’ . 1 . Solar Array Model s 

The I'unlels encountered in the literature are discussed below. Table 
;’-d summarises model capabilities with respect to many attributes, and 
Table I'-^l describes the exist inq modi'l uncertainty/inadequacies. 

The classic first order model for solar cell behavior has received 
much attention and many papers have been published in this a*, ea. Shown in 
fiqure ;’-l is the bC solar cell model. Hie current source 1^^, is an illumi- 
nation dependent constant current soui'ce. Diode D is an ideal diode whose 
characteristics are derived from the 1-V characterist ics of Ihe solar cell 
beinq modeled. Kesistance aiul It^j^ .ire lumped models for various effects 
occurriiuj within the solar cell. If,, accounts primarily for the ohmic resis- 
tance of the sem i conduc t I'l' material of Hie i.'ll as well as the ohmic cont.ict 
resistance ot the cell i:onnei. t ioir. . If, .uiounls primarily lor the surtace 
leakaqe curri'iits around the si'iiiii oiuliic lor Junct ion. 
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MODELS FOR \,7 (RESISTANCE) • ACCEPTADU ACCURACY EOR EARTH ORBIT/ AND SOME INTERPLANETARY, AT 

LOW INTENSITY/ LOW UMP/ CELL TO CELL VARIATIONS ARE NOT PREDIC- 
TABLE. (-.IsC/ '100 C). 


MODELS FOR 2/3 

MODELS FOR <1 

MODELS FOR S 
MODELS FOR 6/9 

MODELS FOR 8 
MODELS FOR 10 

Models for 11 


SMALL SIGNAL AC MODELS AVAILABLE/ NO LARGE SIGNAL OR TRANSIENT 
MODELS AVAILABLE. 

ACCURATE MODELS AVAILABLE. - MAY REQUIRE INTEGRATION WITH PER- 
roRMANCt MODLL. 

ACCEPTABLE OVERALL MODELS. 

ACCEPTABLE MODELS FOR SILICON PIANAR ARRAYS, HOWEVER THERMAL 
CYCl.ING LIT LOTS FOR MISSIONS LONGER THAN '> YEARS IS BASED ON 
EXTRAPOLATIONS FROM SHORTER TERM DATA. 

NO ACCEPTABLE MODELS AVAILABIE. 

SHADOW ANALYSIS CAN BE TIME CONSUMING FOR SOME MISSIONS. AN 
AUTOMAITD SYSIIM IS RFOUIRID TO REDUCE. COST AND ANALYSIS. 

ACCEPTABLE MODELS EXIST. 


CONCENTRATOR ARRAYS 


MODELS ARE BEING DLVILOPLD, TESTING RLQUIRID FOR VALIDATION. 

NILD OPTICAl MODLL 

MIRROR AGING MODEl REQUIRED 


TAHI I 


LOMPONtNl MOlilt EAI'AIUI 11 Y SUMMARY 
SOI AR ARRAYS 


NASA 
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Soi AK AuKAY MuULL 
llNaKTAlNIY/lNAUL«UAi:Y 


• FEW DOCIMENTS DEAL WITH INTEGRATION OF SOLAR ARRAY MODEL WITH ENTIRE 
POWER SYSTEMS MODEL 

• LARGE SIGNAL MODEL REQUIRED 

-- ACCOUNT FOR CAPACITIVE AND INDUCTIVE EFFECTS 
— ACCOUNT FOR RESISTIVE EFFECTS^ CABLING 
— TRANSIENT ANALYSIS 

• PUSMA INTERACTION 

• CONCEiNTKATOK ARRAY OPTICAL MODEL NOT AVAILABLE 


Table 
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Figure 2-1. Classic DC Solar Cell Model 

The advantages of the model pictured in Figure 2-1 are that it is a 
standard model used by many researchers to provide very good first order 
predictions of the DC behavior of the solar cell. The simplicity of the 
model makes it easy to use. It is necessary to obtain only a few circuit 
parameters in order to use the model. Since it is in the form of a circuit 
model comprised of standard circuit components, it can be easily used in 
conjunction with available circuit analysis computer programs. 

Other DC circuit models for the solar cells have been developed. 

The models are generally refinements of the model seen in Figure 2-1, and 
are used to model what might be termed the second order effects in a 
solar cell. Generally, additional circuit elements are added in order to 
produce the d.-sired second order effects. Shown in Figure 2-2 is a multiple 
element solar cell model. This type of multiple element model is used for 
such things as precisely modeling the effects of distributed series resistance 
or complex voltage dependencies. Multiple element models are also used to 
model high illumination (concentrator) devices as well as being used for 

modeling the effects of temperature changes in the various parts of the solar 
cel 1 . 


A ^ 


SC 


Illumination 
dependent, 
const, current 


i 
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Figure 2-2. Multiple-element DC Solar Cell Model. 

The standard large signal AC solar cell model is shown in Figure 2-3. Like 
the DC model, this model is simple and easy to use, producing good first 
order results. Various techniques are given in the literature for measuring 
the data and performing the calculations for obtaining the values of the 
model parameters. For the first order AC model, the parameters are dependent 
upon output voltage and current, cell temperature, cell illumination, and 
of course, physical cell dimension. It is therefore necessary to take these 
dependencies into account when model parameters are calculated. 

I's 

• if — - * -VW. r- 1 




^ — — — 

Constant current source 
Diffusion Capacitance 
Junction Resistance 
Transit ion Capacitance 


R “ Shunt Resistance 
Rj. - Series Resistance 
R| - Load Resistance 
C, Load Capacilattce 


Figure 2-3. Siii.il 1 -s igna 1 AC Sol.ir Ct'll Model 
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Whi'n tlodlhitj with a 'inlar iIIT.iv iimilt'l , «t'. opjur.i'il In Itio snltir ci'l I 
imulol, !.ov(Mm 1 .uidHioiMl cdfi is iiiusi lu' .u i otml oil (im’. Ihosi' oMoi l s .100 
host im liuloil in Itioiiuiilol . 0 . .ultliilonnl i Irniil ('loinonls. Ilii' ml oi’onimoo- 
lion wlrini) ol tlio sol.u’ .u'roy 1*. tmuloh'd h.y .uidil htii.il ro'.i'. loi”. , indinioi’s, 
(Uul iM|Moitoo‘< in both I ho IH' .ind At’ iiiodol'.. ( r.ap.u’ i I iim o and indui tani o 
havo no oiloot in llio lUl inodol.) Iho valui's ioo I hoso aildilional iiiodol 
paraniohM’s ai’O ilopondoni on an’ay •. i/o. 

Suhsh'ato parasilio i apao 1 1 ami' oan hoooim> n iiini i'ioanl and Ihoroforo 
must ho acoountod too in llu’ .olao aooa.y AC modi' 1 . Induotaiuo ol'toois poo- 
ducod h.y Iho inoidont illumination ( phot oinduo tamo) liavo oocontl.y boon 
slu)wn to ho on tho oodoo of 100 pll/om:’ wind) oan ho significant in a solar 
aooa.y. llio valuo ol plioioiiului lanoo olisoovod is a function of I'oth illumina- 
tion intonsit.y and illumination wavoloiuith. Mooo work noods to ho dono to 
pooduco an appoopoiato cioouit olomont for inclusion of this offoct. 

In oodoo to poodicl tito d.ynamic offocts of cliaiuios in tho solar aooa.y 
such as coll failuoo, shadinii, oclipsinq, 00 load cliaiuios, it is nocossao.y 
to uso a laogo-siiinal modol. A satisfactory laoi|o-signal model for llio solar 
coll has not boon found, i’oovious dynamic models usod tlio paoamotoos C^ , 
t'd* Rj. R^, and from tlio AO modol to form a transient modol. However, 
hocauso of tho siijnificant variations in those paoamotoos over wide oatuios 
of voltaqos, currents, illumination and tompooatuoos, this oxtonsion of 
tho AC model is not suitable as a lariie-siqnal model. A moans must ho 
found to incorporate tho appropriate circuit effects over the full large- 
siiinal ratuio. Tho typical lumped parameter modol iiui technique may prove 
to ho difficult to uso lor this task. 

;M . . 1 Array Model i ng Approa ches 

Tho literature survey determined that there is a large amount of 
material dealinq with the various aspects ol solar colls, however, few 
documents deal with the integration of tho solar coll model with the power 
subsystem. 
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Several approaches are proposed for modeling a solar array, first, 

It Is possible to model an array by an Interconnection of individual solar 
cells (individual cell model approach) ns shown In Figure 2-4. This tech- 
nique Is straightforward and does not require a lot, of model Ing experience 
1n order to use it. The parameters for the Individual cell model are 
readily obtained using established measureinent/calculation procedures. 
Shadowing effects as well as cell faults can easily be Inserted Into the 
array using this modeling approach. Since the individual parameters for 
each of the cells is entered into the model, the effects of parameter 
variation over an array of cells can readily be included. Additionally, 
changes in array structure as well as array expansions can be made by 
simply specifying a different interconnection scheme. Thus, the individual 
cell model approach seems to bo quite general and conceptually easy to use. 

The individual cell model does have certain disadvantages especially 
when larger arrays are treated. Most circuit analysis programs are not 
capable of analyzing arrays larger than a hundred cells. Those that will 
handle the larger arrays use an extremely large amount of computation time 
in performing the analysis. In addition, the previously mentioned advantage 
of being able to enter and vary each cell parameter on an Individual basis 
can oe a disadvantage In terms of user time for large arrays. 



Figure 2-4. Individual cell approach to modeling a solar array. 
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Anothor rtpnro»>t'h is to os(' o sint]li' loll iiuulol to sinuilatc' tho ontiro 
array (macro inodol approach). Soi' ri<|nr(' iho major advantaqo to this 

approach Is lis ahillty to m1ii1mi/.(' coiiipoli'r Mine aod space oecchsl to 
analyse array hidiavior, ParamoLers can tie ohtaiiK'd In sevc’ral way. from 
Incurporatlnp tho moasiirement/calciilatlon of l hi' individual cells accordinti 
to tho Inti'rconnecI ioir. to inakim.) lerininai iiiea'.nremenls <it tho array'f^ 
connoctlon points. Neither of these nii'lhoils lakes much more effort than 
would he needed for the individual cell model. Once these paramelors are 
determined, only a simjle set of data is needl'd for inpul to i he analysis 
pro 11 ram. 

The disadvantages to tho macro model approach Include the fact that 
a new set of macro model parameters needs to he developed for each different 
array configuration. This tends to make the model loss flexihlo than is 
desirable. Furthermore, shading and individual cell faults cannot be 
easily included in the macro model without, again, a redevelopment of 
model parameters. Individual cell parameter variation studies are also 
difficult to perform using the macro model approach. 



} 


Figure 2-5. Macro model approach to modeling a solar array. 

An approach which uses a combination of the macro model and the 
individual cell model as shown in Figure ?-6 could be very advantayecus. 
This approach retains the flexibility of the individual cell model approach 
while still reducing tho necessary computer time and space requirements. 
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of 



Coll fault..-, and groupri of coll fault.r, can bo r,)inulatocl with tho combination 
model, Shadowl.ig offoctr, can bo accounted for in much tho samo way a> would 
bo dono In tho Individual niodol . Pnrainotor variation '.tudles can alr.o bo 
performod by varying Individual coll parainotorf, of I ho toniblna- 
tion niodol and leaving tho macro modol nuhnoctlon paraiiiotorr, fixed. 

Tho disadvanl.ago of using tho combination modo! In tho need to 
doterinine at Ic'ant two sots of modol paraniotors-»ono not for the group of 
Individual colls and ono sot for tho group of macro models. Another 
disadvantage Is the need to dovolop a now sot of macro parameters when It 
is desirable to Include more or loss of the Individual cells of the array. 


icell 

1 

cell ! 

1 

r- 

i cell 

1 

|cell 

i 

1 

1 



r 




1 r 

1 . . 

cell 


macro model -subsection 


Figure 2-6. Combination macro model and individual cell model. 

The choice of the modeling approach used depends upon the application 
at hand. For example, the macro model approach would be used for performing 
an overall system analysis. In this type of analysis it is presumed that 
the model of the solar array need not be as flexible or general as other 
types of analyses. For an analysis which is to investigate system component 
details, e.g., battery conditioner, etc., the solar array macro model will 
suffice. The individual cell combination model would be used when Investi- 
gating solar array details. For this typo of analysis It is likely that 
macro modol s would be used for the other system uimpotioiits , 
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?.],?. ,Z So] u r A r r ft y Mod o 1 R n t^i i ij ’_oiii on t r. 
Inip ortft n t Pft rft iiin torn 


Tho ob.inctivn n1 tlilr, nnctlon in to dntorniino a proriiiilriftry not of 
porformanco ntandartlr. and oxpoctatioiis wlridi tho nolar array inodol niunt 
satisfy. Important solar array model output paramotors Include tho array 
voltacjo and current and array power delivered. Temperature variation 
Information Is a necessary model output parameter. Temperature Information 
needs to be coupled to tlie rest of tho system to determine system heat 
flow. A more comprehensive model iiuj package would allow tho user to 
perform further analyses on solar array size, magnetic effects, radiation 
effects, plasma effects, as well as aging and other degradation effects. 


The required Inputs for the solar array model Include a set of 
parameters for each of the elements of the array model being used. When 
using either the combination model or the individual cell model. It is 
necessary to input the interconnection scheme for the array. Another 
important input parameter is that of the array illumination characteris- 
tics. Included In the illuminatior. characteristics are such items as 
intensity, frequency, and shadowing information. The ambient temperature 
must be included as a model input parameter in order to provide an initial 
point for a thermal effects analysis. For analysis which will perform a 
study of the behavior of the array when subjected to fault conditions. 

a set of cell faults must be specified. These faults should be of both 
the fixed and time- varying types. 

The way in which the input parameters are entered depends upon the 
nature of the parameters as well as the way in which the parameter is 
measured/calculated. Parameters which are single valued can be input as 
con..tants. Other parameters which may be functions of network variables 
can he input in two different ways. Functional descriptions can be used 
to specify parameters whose values fit into particular functional forms. 
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fl •.,>co„d W.IV In wlUnli this ty|.n of |.„,-„i«o(,t cm ho Input |,; „„ , ,(.,1 of 

data points. TIic'm. data points am llkoiy to iiavo |,ono nlitalood from array 

«asur„,,onts. Valuos nsod l,ni„™n data points cai, oltluT ho oonstant at 

tin last spoclflod data point vain,, or ..in ho a vaino ohtalnod by intorpo- 
Ifitlori. ' ■ 

J OillvJ 1 1 ,1 onion t. R, 0 (,|u 1 r onion 1, s 

In ordor to no,.roct1y mod,, I tho solar ari-ay, U,o liaidollnp roqulro- 
nionts of tho Individual oloiiionts inio.t |,o oxamlnod. tt H essential to use 
the cori-ect 1-V chai'acterfstlcs for all clonients of the model. Importantly 
t ementb which arc dependent upon other circuit parameters must exhibit tho 
proper dependency cha,.acter1sttcs. Those characteristics must be determined 
n the measureiiient/calculatlon phase of the model development, lor large 
a, -rays. It is very Important to Include the subst,-ate capacitance effects. 
Careful examination of the resistance, Inductance, and capacitance of the 
Interconnections structure is required. In addition. It is neccssar.,. that 
each Of the individual parts of the solar arrav model be capable of pro- 
ducinq infornidtion usable in a thermal effects analysis. 

2 . 1 , 2 . 3 Model Development Met hn d ole gy 

A solar array model development should start by first devising an 
appropriate circuit model and element configuration. During this phase of 
he development. It is necessary to determine a set of requ-red parameters 
frr a DC. AC, and ti-ansient analysis of an array. The parameter set for 
K Individual cell model would be developed first, followed by the macro 
model and combination model parameter sets. Included In this phase of 

the development Is the incorporation of the ehmient parameter tnterdepen- 

cJencies, 

The next phase of the development involves the verification of the 
model using experimental data. Model responses for typical and worst-case 
-X "f'Pd to be calculated and compared to measured data, limitations 

on the bounds of worst-case situalions are determined by the experimental 
capabilities of the model developers, broadening of tho worst-case bounds 
as much as possible allows for more confidence in model predictions, 
produced for situations outside the measurement bound. 
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(1110 (111(1 two IldVO lioon MK-COS'.I 111 1 V Mtlll|l] 0l0(| t lio iiiodol 
(’ill) ho iifiod III priHlIct r,o|,ir iiriMy hohnvinr. II Is n\ till:, polnl. I;ha1; tho 
ittodol is inr.nrpnr/it.oil into thi^ ovoral l powi'r s.ystont niodol . 

?.A ,?A M.a t homa 1 1 ca 1 Skills l<o()ii i r od 

Tlio (lovolopoi* of till) solar array inodol imisl. havo a bafric. oloctrltal 
engineer lug knowledge of spacecraft power generation subsysti’ins and in 
particular, a significant knowludije of solar cells. Furthoriiiore, the de- 
veloper must have an understanding of circuit element representation and 
of the interdependencies of physical effects, e.tj., the dependency of Cj 
on junction voltage. This requires a familiarity witii semiconductor 
electronics. 

In order to aptly implement the solar array model and integrate it 
within the comprehensive power system analysis package, it is necessary for 
the developer to have some basic programmer skills, and some understanding 
of numerical analysis methods. A background knowledge of available circuit 
analysis programs such as SPICE or SCEPTRE would also be very useful. 

2 . 1 . 2 . 5 Modeling Objectives and Adequacy 

The prime objectives of the solar array model is to accurately 
predict the behavior of an array of photovoltaic cells, and be usable 
in a comprehensive power system model. 

The solar array model must be capable of performing in three different 
analysis inodes: DC analysis, AC analysis, and transient analysis. The 

individual objectives in each of the three different analysis modes are 
listed below: 

DC “ voltage vs. illumination 

- current vs. illumination 

- output power vs. illumination 

- teniporature cffecis predictions 

- load variation prediction 

- sensitivity etfects due to resistance changes 

“ fault condition predic. l iinis 
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AC - oiitfuil. (wr’diiH’tf'r''. v^. . Trni|Ui n( y 
-- outRUt dol.ni’iirliKn: Ion 

r.onr.lt ivUy nffnctr, diir' to R’s , L’f. , (iiid C'n 

" fdiddow’lno 

- foult: trAnr.Inntii 

- triMvijonl!, 

A &ot of nfltoriu are dufinod liulow wIim.Ii can Uu used to dctcnnlno 
tho degrof! of adequacy of the liolar array iiiodol : 

• Accounts, for all physical effects including the proper I-V 
characteristics, Impedance, sluidowlng effects, fault effects, 
in all applicable modes. 

• Obtains the necessary degree of precision, 

• Performs the analysis necessary without using Inordinate amounts 
of computer time or space. 

• Usable In conjunction with other power system component models. 

• Does not require a high degree of mathematical skill of complex 
ineasurements/calculation so that it is easily usable. 

2 , 1 . 2 , 6 Necessary Developmen t_ £' f f o r t^ 

Two areas of development effort are necessary in order to produce a 
solar array model which can be used in the comprehensive power system model. 
First, the presently available model needs to be improved to satisfy the 
previously stated objectives. 

The development efforts which still need to be performed fot 
the presently available mode! (individual cell model) are; 

• Account for all series resistance effects within the array. 

• Account for all inductive and capacitance effects within the 
array. 

• A means of methodically specifying various interconnect schemes. 

• A means of easily accounting for the effects cf shadowi.ig and 
cell faults. 
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The second effort necessary is to develop an appropriate macro model 
which satisfies the objectives. Third, a method must bo developed for de- 
termining the parameters for the many possible types of macro model subsections 
which will bo used in the combination model. Besides development of the macro 
model subsection itself, the other important components of the model which 
must be developed are: 

• The interface between the macro model which will convey other 
important array information, e.g,, thermal effects Information. 

• An accurate set of parameters for macro models for various 
sized array subsections. 

• A method to input shadowing and fault data. 

• A means of analyzing the sensitivity of the array performance 
to parameter variations. 

In addition to the development efforts listed above, pertaining 
primarily to the solar array, it is necessary to develop a structure for the 
overall system model into which the array model will fit. This structure 
may be loosely defined initially, but will likely become more constrained 
as the development of the comprehensive model system progresses. 

2 . 1 . 2 . 7 Uncertainties 

The major uncertainty for the solar array modeling effort lies in 
the transient analysis model. Little is reported in the literature on a 
working large-signal model for solar arrays. The amount of time and effort 
necessary to extend the DC and AC models to a workable large-signal model 
could be substantial . 

The otlier area of uncertainty in the solar model deyelopment is 
that of the macro model. The points which seem most uncertain about this 
development effort are: choice of the appropriate parameters, method for 

calculating/measuring macro model parameters, and how can flexibility be 
built into the macro model or combination model scheme. 
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.3 Shunt, Rogulotor 

Thorc wuro fow pub 1 leal ions oti shunt roqulator niodols. Ilio Utora- 
turo usually addrossos tochnigut’s and design ot shunt eegulators, and not 
coniputor modeling sr.heinos. Tal)le J’-b suniinarizos the present capabilities 
of shunt regulator models with respect to various attributes. Table T’-6 
summarizes the uncertainties/inadeguacies of existing models. 


Schematic and Definition 


*R-, Bus Voltage '^b 
o \'V\ —f 


«1 <Vs 



INPUTS 

Typically 

OUT PUTS 

Typically v^. ^ i^, & 


*The term shunt regulator inherenlly includes a st*ries droppittg resistor, R^. 
Since this resistaiut’ is the -.olar array series resis.tance in this case, 
circuitry excluding thii. element is I ( rmed a shunt limiter sime it only 
limits bus voltage upper value. 


TWO CA TEGORIES 

• Linear 

• PWM 



SHUNT REGULATOR 
MODEL ATTRIEWTES 

MODELS FOR 1>2 

MODELS FOR 3,5>6 

MODELS FOR 
ALL MODELS 

Table 2-5. 
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WELL UNDERSTOOD AND ACCURATE FOR SINGLE ARRAY LUMPED SHUNT 
NO MODELS FOR SASU AND PARTIAL SHUNT * 

LIMITED CAPABILITY FOR SINGLE ARRAY LUMPED SHUNT 
NO CAPABILITY FOR SOLAR ArRAY SWITCHING UNIT (SASU). 

NO MODELS AVAILABLE 

INHERENTLY LINKED TO SOLAR ARRAY. DECOUPLING IS EASY 
FOR SINGLE ARRAY SHUNT.. BUT MUCH MORE COMPLEX FOR SASU 
OR PARTIAL SHUNT CONFIGURATION, 


COMPONENT MODEL CAPABILITY SUMMARY 
SHUNT REGULATORS 


) 


2(i 




NASA CR" 166820 
TRW 38651-0001 


Shunt Regulator Models 
Uncertainty/Inadequacy 


• VERY LIMITED PUBLICATIONS ON SHUNT REGUUTOR MODELS 

- DESIGN ORIENTED 

- INCUIDED IN SOME OF THE OVERALL MODELS^ PSIM, ETC, 

• GENERALIZED CIRCUIT «^'r‘EL CAN BE UTILIZED 

~ SPICE HAS BEEN USED BY TRW FOR ANALYSIS OF SHUNT REGULATOR/LOAD 
filter interaction AND STABILITY MARGIN ANALYSIS, 

- verification BY BREADBOARD AND ENGINEERING MODEL TESTS (UNPUBLISHED) 

- SIGNIFICANT ANALYTICAL EFFORT REQUIRED TO ESTABLISH THE MODEL 

- NO POWER DISSIPATION OR ELEMENT FAILURE ANALYSIS 

• NO DOCUMENTED TRANSIENT MODELS WERE FOUND 


Table 2-6 
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PMlOJ t . Oporati ona l Modo;; 

Mode 1; Sensed voUayo below setpoint voltayo to activate shunt. 

Shunt active olemopt off. 

Mode i!: Sensed voltage above setpoint voltayo. Shunt elements 

activated. Linear type shunt lias some active elements 
in linear region. PWM shunt at less than maximum duty 
cycl e. 

Mode 3: Sensed voltage above setpoint voltage. All linear shunt 

active elements saturated, PWM shunt at maximum duty 
cycle. 

p^tprmance 

The standard criteria for evaluating the performance of the shunt 
regulator are as follows: 

• DC Voltage Regulation* 

• Transient Response* 

• Stability 

• Starting 

• Saturation 

• Active Element Dissipation 

• Output Impedance (dVj^/dig) 

• Regulation Sensitivity (Vj^ max - V|^ min) 

• Input Regulation* x V|^)) x lOO';- 

• Load Regulation* min load - full load)/Vb min load) x 100 

• Temperature Coefficient*, TCC:./C") =• ± (V,, max - \ min) x lOOW 

• Ripple Rejection* 

^ Bus voltage 
1 Shunt ( iiffc'Mt 

V - Solar array open ^ iruiit voltage 
0 ( 

*Performance evaluation must consider solar array open circuit voltage 
and series resistance (V,^ fti R, ,). 

()L • > 1 V * 
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? . 1 . 3 . 1 Model ing A pproaches 

Two modeling approaches have been analyzed, 
approach and the functional block approach, 
are shown below. 


Network Approach 

1. System expressed 1n terms of 
Individual elements. Each 
element usually defined In 
terms of a thru variable & 
an across variable. 

2. Basic elements are storage ele- 
ments, dissipative elements, 
sources. Allows use of standard- 
ized models. 

3. Good component level visualiza- 
tion. 

4. Requires high level program, large 
computer capacity, & numerical 
analysis skills In operator. Pro- 
gram usually performs as network 
analyst liiternally. 

5. Easy Initialization operations. 

6. Easy to alter parameters of net- 
work for sensitivity analysis & 
worst case analysis & design change. 

7. Network may be very large at 
system level. 

8. Some programs may require further 
data reduction to obtain frequency 
domain data. 

9. May be able to handle multiple 
topology & discontinuous systems. 


The_ are the network 
The characteristics of each 


Functional Block App roach 

1. System expressed In terms of block 
diagram, differential equations, 
or transfer functions. 


2. Basic elements are mathematical & 
logical operations. System models 
frequently custom developed. 


3. Good system level visualization. 


4. Requires network analyst to 
develop models as functional 
blocks. 


5. Substantial initialization effort. 

6. Difficult to alter parameters as 
required by sensitivity analysis 
or design change. 

7. Hard to handle discontinuous systems. 


8. Must provide all Initial conditions 
for state variables. 


9. No need for equivalent circuit. 

10. Most useful In frequency domain & 
small signal transient analysis. 
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MATHtMATICAL AND KNOlNhtRlMG SKILLS REQUIRLO TO DCVCLOP 

SHUNT REGULATOR MODEL 

Tho dovolopniont of .1 shunt i'oi|ulntor niodol roqiilros basic oloctrical 
ongineoring knowlodgo of spacocraft power gntioral ion and distribution 
subsystems. In addition, tho following skills are required: 

• basic electrical engineering knowledge of spacocraft power 
g e n e r a t i 0 ti s u b s y s tom . 

• Network analysis skills. 

• Semiconductor modeling and analysis skills. 

• Numerical analysis and computer programming skills. 

• Non-linear analysis skills. 

Important Parameters to be Analyzed 
1 . DC Model ing 

• Worst-case: drift, eclipse, radiation degradation, mission life 

fatigue, temperature variation, etc. 

• Voltage regulation under all variations of bus loads and solar 
array configurations. 

• Survivability: ability to continue regulation in the event of 

component failure. 

t Semiconductor component stress. 

• Shunt contribution to spacecraft thermal loads. 

• Power regulation of peak power tracker (PPT). 

2* Small S ign al AC Model 

• Stability under worst-case as listed in DC modeling. Quantify 
stability in terms of gain and phase margins in the frequency 
domain against worst-case variants. 

• Output impedance versus frequency. 

• Interaction with system components: bus filter, switching regu- 

lator input filter, solar array switching, switched shunts. 

f Shunt element overlap: increasing sensed voltage causing addi- 

tional shunt element to become active and decreasing sensed 
voltage causing saturated shunt element to come out of saturation. 
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§ Small ampllt.udo stop or »’amp tranr.iont rosponso. 
t Iruilvldual control loop oluiractorlstlcs. 

• Operation at threshold of solar array section switching. 

• Stability of variable reference voltage lontrol loop for peak 
power tracker operation. 

« Phase relation between sensed voltage (regulated voltage) and 
shunt element voltage for partial shunt. 

3. Larg e Signal Mo del 

The large signal model must accurately predict performance during 
step load changes, solar array section switching, surge suppression 
(such as solar array radiation barrage), and setpoint voltage cross- 
ing. 

The effects of the following parameters must be considered in a 
transient response model. 

• Regulation sf nsitivity. 

• Sequential shunt switching. 

• EMI response 

• Non-linear characteristics. 

• Settling time. 

• Sensed (regulated) voltage overshoot. 

• Sequential element burn-out. 

t For PWM shunt: steady-state ripple injected onto the bus. 

^ • System I nt erac ti on 

The shunt model must account for the following interactions with 
the power system. 

• Interaction between shunt (bus) filter capacitance and resistance 
and the shunt. 

• Interaction between switching regulator input filter and the shunt. 

• Effect of step load and source changes (within and beyond the 
limits of shunt sens ; ; ivity) . 

• Effect of periodic load and source variations both sinusoidal 
and non-sinusoidal (within and beyond the limits of shunt sensi- 
tivity. 

• Interaction ot stiunt with mniative resistive characteristic ot 
the power conditioning switching regulator. 

• Effect of shunt on solar array operating in current /vo 1 tage model. 

:'-3i 
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MODELING REQUIREMENTS FOR INDIVIDUAL PIECE PARTS USED IN 
SHUNT REGULATOR MODEL 

In some cases a transfer function approach Is not adequate to predict 
overall performance, and It is necessary to model some of the parts or 
functions internal to the shunt. Modeling requirements are given as a 
function of the shunt mode, 

1 . Mode 1 Operatio n 

Ne twork Appro ach. Shunt filter element values, bus loads, solar 
array segment mo'dels. 

Functional Block Approach , Shunt filter transfer function, solar 
array source mathematical description, and load Input Impedance 
description. 

2. Mode 2 Operation 

Network Approach . Solar array source: Topic discussed under solar 

array modeling and analysis. 


Reference voltage: Independent or dependent voltage source equiva- 

lent to actual zener or peak power transfer characteristic. 

Error amplifiers: Feedback circuit passive elements values and 

amp model element parameters including input capacitance and 
resistance, DC gain output voltage dependent resistance divider 
expressions. 

Majority voting: Generalized model to be developed. 

Shunt element seauencer: Generalized model to be developed. 

Sequential shunts: passive component values, parameter evaluation 

procedure for voltage dependent capacitors and current sources of 
active shunt elements. 

Bus filter: passive element values capacitance, capacitor equiv- 

alent series resistance, swamping resistance. 

Load model : As determined by applicable modeling and analysis 

section dealing with particular load. 

Functional Block Approach . 

Transfer functions: Determine range of linearity about appropriate 

operating points for AC small signal models for each functional block 
as previously listed. Determine minimal mathematical expression to 
predict gain and phase characteristics of each functional block. 
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nifferentlal equations: Determine linear and non-linear differen- 

tial equation for each functional block. 

Mode 3 Qpor^lon 

Shunt filter element values, bus loads, solar array 
segment models, and shunt element dissipation resistor values, 

BlM|_]pnaX.j^l(K:k Transfer function of shunt filter In 

parallel with shunt element dissipative resistors, solar array source 
description, and load Input description. ^ 

MODELING OBJECTIVES 

1. Test performance of shunt at system and subsystem level against 
design requirements. ^ 

tioJlI^^ sufficient flexibility to accept various shunt configure- 

3. Serve as design aid to analyze existing designs and help effect 
Unproved design. 

power generation subsystem components mission performance 
1n advance of hardware construction. 

5. Util Ip existing models and analysis programs to the maximum 
extent possible. 

MODELING ADEQUACY 

1 . Network Approach 

• Modeling techniques for components exist. 

• Computer programs exist. 

• Can provide analysis objectives of AC and DC modeling. 

• Easily aaapted to new technology. 

• Exact models not yet developed. 

• Cataloged procedures for modeling not yet developed. 

!. Functional Block Approach 

• imited existing modeling techniques, 

• Models not easily adapted to new technology. 

• pnoot provide verification of design criteria such as power 
dissipation^ component degradation, etc. 
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• Coinputor anil lysis programs ox1st, 

• Dev.^lopmont efforts to satisfy modeling objectives are quite 
imposing and must be done individually for each of the typical 
shunt regulator configurations. 
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I 

^ • 1 • A J^P.nj/.e r t 

MoTit of t.ho (OC-DC Convortorr,) iiindnlr, clnr,crll)od In thir. ".ocflon woro 
dovolopod to analyzo or r.1niulal:o a voltatjo rnqulator from l.ho point of vlow 
of convortor porfortnanco, and convortor oporatlon (modo). Stability analy" ■ 

sis was the prime motivation for most of tho work reported In the literature. 

As models were developed, by-products became available, such as audio ; 

susceptibility, output impedance, etc, However, the kind of model that Is 

il 

needed for power system modeling is completely different than the ones j 

developed to date. j 

Table 2-7 summarizes the capabilities of present DC-DC converter i 

modeling with respect to various attributes. Table 2-8 summarizes the j 

existing models uncertainty/inadequacy. 

The linear regulator has the advantage of ease of design and analysis, > 

and no EMI problems. The disadvantage is that it is dissipative and heavy. j 

The switching converter has less mass and less thermal dissipation than 
the linear regulator. It is more difficult to design and produces EMI. j 

2 . 1 . 4 . 1 Switching Regulator Performance Categories 

1. DC Model Performance Categories: 

input voltage range , 

output voltage range 
output power range 

DC regulation , 

mode? of operation - continuous current, discontinuous current 
semiconductor component stress 

I 

worst-case analysis i 

EMI performance j 

output ripple voltage | 

2. Small Signal AC Model Performance Categories: | 

stability (margins of stability) 
audiosusceptibil ity 

output impedance i 

step load transient response (small load transient) ■ 

i 

input filter interaction i 
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single and multi- loop control 

duty cycio pul no modulation nchomor. 

charactorlstlcG and componGatlon of control loops 


3. I.argo Signal Modol Porforinanco Catogorlcv, 

stop chango of Input voltago 
stop chango of oiitpiii currniit 
global stability (largo ".liinal stability) 
starting transient 

noni inearltles such as saturation of magnetic components, 
saturation of Op Amp 

protection circuits 

4. System Interaction 

Interaction between switching regulator and Input filter 

Interaction between switching regulator and unknown source 
Impedance 

interaction between switching regulator and payload 
(reactive load) 

2 . 1 . 4 . 2 Switching Regulator Modeling Techniques 
Average model (Caltech ) 


Power Stage: Taking advantage of the much lower output filter 

resonant frequency In relation to the converter 
switching frequency, the nonlinear switching power 
stage is approximated by a continuous small signal 
1 Inear model . 

Approaches : 

• topology deduction to form a linear circuit model 

• equation derivation to form a linear state-space 
model 


Both applicable to either continuous current opera- 
tion or discontinuous current operation 

Analog to Obtain output-voltage-to-duty cycle transfer function 

Discrete-time through describing function technique. 

Conversion : 


Input Filter: Identity interaction between the input filter and 

other two tunctional blocks. 
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Component Model Capability Giimmarv 
nC"DC Converter 



MODELS POR 1,2 
MODELS POR 3,S 

MODELS FOR 
MODELS FOR 6 


WELL UNDERSTOOD^ ACCURATE 
LIMITED CAPABILITY 

- PRECISE INPUT DATA IS REQUIRED POR LARGE SIGNAL AC ANALYSIS 
LIMITED CAPABILITY 

WORST-CASE MODELS EXIST 

- PRECISE INPUT DATA REQUIRED 


Table ?-7 
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DC-DC Convert tiR 
Uncertainty/ Inadequacy 

• DC - S^«LL SIGNAL AC MODEL 

• DISCRETE AVERAGE MCX)EL 

■ WTA BA% DOES fWT CONTAIN ALL COMMONLY USED CONFIGURATIONS 

- CURRENTLY USED VALIDATION TECrtllQUE LACKS ACCURACY 

• DISCRETE TIME DOMAIN MODEL 

- COMPLEX, HIGH COST 

- COMPONENT LIBRA'^Y NEEDS EXPANSION 

• LARGE SIGNAL MODEL 

• AVERAGE M(X)EL 

- LACKS ACCURACY ON SHORT TIIC CONSTANT TRANSIENT 

■ PROPERTIES OF OP-AMP PROTECTION CIRCUIT AND TRANSFORMER 

wMIUKAl Iv^ 

• DISCRETE TIME DOMAIN MODEL 

• DIFFICULT TO GENERALIZE FOR UNIVERSAL APPLICATION 

• GENERAL CIRCUIT ANALYSIS PROGRAMS 

• DETAILED ANALYSIS RESULTS IN CUMBERSOME INPUT AND HIGH RUN COSTS 


Table 2-8 
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OF pOOn (.RlAMi'V 


Merits 

The average model Is simple and easy to use, and It Is readily applicable 
to complex circuits and systems. The model addresses both DC performance 
and small signal AC low frequency performance. It accepts either the 
transfer function form or the equivalent circuit form for the linear model. 

Limitations 

The average model only preserves the input-output properties of the con- 
verter. The original properties of a state variable (Inductor current) 
is lost In the process of averaging. Consequently, the average model 
cannot be directly used for multiloop control systems which sense the 
output voltage and Inductor current (or voltage). 

The large signal model is not readily available for transient and start- 
up analysis. There 1s diminishing accuracy beyond 10-20% of switching 
frequency. Not suitable for high-gain wide bandwidth regulators, such as 
high performance regulators employing multi -loop control schemes. 


The canonical circuit model cannot be used directly to implement the 
multi-loop control. 

Exte nded Average Model 

The extended average model was developed for application to multiloop 
control systems. The three-input, two-output average model Is shown below. 



SMALL SIGNAL POWER STAGE MODEL 



POWER 

STAGE 

MODEL 




Inputs: line disturbance v^ 

load disturbance 1^ 

duty cycle disturbance d 


Outputs: output voltage disturbance 
switch current disturbance 

P 


2-39 


NASA CR-1 66830 
TRW 38651-0001 


Tho QxtsndBd niodol has the following added advantages ovor the 
original average model ; 

1. Equivalent circuit model representing the original properties 
of the converter's Input, output, state variables. 

2. State variable model for state space analysis. 

3. Transfer function model for classical frequency domain 
feedback control analysis and design. 

Discrete Im pul se R esponse Moi^l 

This method Is based on calculation of the state vector perturbations 
at the switching Instances; these are discrete in nature. The 
approach Is to represent the power stage with a linearized discrete 
impulse response model, and then applying the z-transformatlon. The 
discrete time domain model is then transferred into the frequency domain. 
Merits 

Accurate power stage model up to one-half the switching frequency for 
AC analysis. (The theoretical limi', of any linearized model.) 

Limitations 

The model requires complex analytical derivations, requiring a high 
degree of mathematical background. It is difficult to incorporate an 
Input filter, and the model only provides duty cycle to output voltage 
transfer function. The model cannot be readily used to study disturbance 
from the line voltage and the load. 

For a complex switching converter system, no closed form analytical 
model can be derived. Numerical techniques have to be comployed. 

Discrete A verage Model 

The discrete average model is derived by combining the techniques of 
the average model and the discrete impulse response model. The approach 
is to use the average techniques to derive a model for state variables. 

(State variables are well behaved and continuous.) Then, use discrete 
time representation to derive an output voltage expression. (Since the 
output voltage is discontinuous due to the filter ESR.) 
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Merits 

Improved accuracy of the average model In high modulation frequency 
up to one half of the switching frequency. The technique retains the 
simplicity of the average model. The model Is easy to use, and Is 
presented In the form of circuit model, transfer function model, and 
state space model. It Is suitable for single-loop and multi-loop 
control modeling and analysis. 

Limitations 

Difficult to use for transient and start-up analysis. 

Discr ete-Ti me Doma1>i Analysis 

Using state space representation, a nonlinear discrete-time system Is 
derived that models the converter exactly. This system Is linearized about 
Its steady-state solution. 

Modeling Approach 

Exact formulation of state equations. Use Newton Iteration to solve 
for the exact equivalent state. 

The system is linearized about the equilibrium to obtain linear 
time-invariant model. 

Use z-transformatlon and obtain frequency domain transfer function 
representation. 

Merits 

No assumption is made. Most accurate small signal linear model for 
stabll Ity analysis. 

Can predict high frequency (subharmonics of the switching frequency) 
instabll Ity. 

The formulated recurrent state equation leads to a cost-effective 
performance analysis. 

Limitations 

Basically, it is a numerical analysis. No closed-form solution Is 
derived to provide physical insight. It is a small signal model, and 
requires background in state space modeling and numerical analysis. 
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Discrete Time Doma i n Simulation 

The discrete tlrue domain analysis techniques are extended to large 
signal simulation. 

Modeling Approach 

Exact formulation of state equation. 

Based on recurrent discrete time domain analytical expression. 
Propagate recurrent equation through numerical computation. 

Merits 

Large signal analysis and simulation. 

Exact duplication of the circuit behavior. 

Capable of start-up simulation. 

Simulation of large signal transients. 

Capable of incorporating all system non! Inearltles. 

A combined analytical and numerical scheme that provides a cost- 
effective simulation faster than other general purpose simulation 
programs such as SCEPTRE, SPICE ICAP, etc. 

Limitations 

Relatively Ineffective for large system simulation, since the user 
has to provide state space representation of the system. 

Large Signal Average Model 

In this model the average techniques previously described are extended 
to perform large signal analysis. 

Power stage model: Using the average technique to represent the 

switching power stage by an averaged continuous-time domain equiva- 
lent circuit. 

Error processor model: Since the error processor Is a linear, no 

approximation Is made. The saturation effect of the op-amp how- 
ever can be incorporated In the model. 

Duty cycle pulse modulator model: The exact duty cycle implementa- 

tion can be simulated incorporating the basic ramp and threshold 
implementation. 
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Merits 

The simplified power stage model provides an effective large signal 
simulation, 

The model is easily adaptable to existing circuit analysis programs 
such as ICAP. SPICE 2, SCEPTRE, etc. 

Limitations 

In general, It Is difficult to Include some protection features such 
as transistor peak current protection, because the inductor current 
(and transistor current) Is approximated by its average value. 

It Is difficult (If not Impossible) to be extended to simulation of 
converter employing a multi-loop control technique where the Instan- 
taneous Inductor current (transistor current) or Inductor voltage Is 
sensed to provide the necessary ramp for duty cycle Implementation. 

It is difficult to include different duty-cycle control schemes such 

as constant T^^ control, constant control, constant frequency 

control, and variable T , variable T .. control, etc. 

on off 

A Combined Discrete and Average Technique for Large Signal Model 
Approach 

A compromise between complexity and accuracy. The converter is first 
represented by discrete time equation. The system is then approximated 
by a continuous time representation that remains in the nonlinear proper- 
ties of the original system. 

Merits 

Nonlinear time varying circuit for large signal simulation. 

Linear time invariant circuit model for small signal analysis. 

Capable of Implementing different duty cycle control. 

Capable of Implementing single-loop control and multi-loop control. 

The model remains in the nonlinear properties of the original system 
and therefore is able to Implement the protection features and 
various saturation effects of the system. 

Liini tations 

The method Is not well documented, therefore the utilities and limita- 
tions of the method are no. clearly understood at the present time. 
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Circuit Analysis and Sim ul ation Programs 

The direct simulation method Is capable of giving detailed Information 
about the system because of the detailed modeling associated with each elec- 
tronic part. 

Approach 

Use existing circuit analysis programs such as SPICE 2, SCEPTRE, ICAP, 
etc., to simulate switching power converters. 

The component models such as semiconductor devices and magnetics are 
presented either In the form of equivalent circuits or state equations. 

Merits 

Easy to Implement provided that the component models are available. 

Easy to determine the component stress and circuit particulars. In 
general. It Is more versatile and flexible to accommodate changes In 
circuit parameter values control modes, and output options for a large 
class of converter circuits. 

SCEPTRE program Is capable of Interfacing with any FORTRAN subprogram 
provided by users. 

Limitations 

Lack of semiconductor and magnetic component library for switching 
applications. 

Lack of effective numerical Integration routine to handle stiff 
differential equations with wide separation of time constant. This 
Is often the case when non-ideal switching component models are 
utilized. It results in either excessive computer execution time 
or numerical instability. 
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6. Cost -I’l Ttu:t i .(’ru’s , 

/. |.ji.y Lu II SI.' 

0. Uonoral ity/1 1niHnt1on to porticiilnr powor stage configuration 
9, Gonoral ity/1 Imitation to particular duty cyclo control schnnio 

10. Gonoral ity/1 imitation to convortnr operating mocio 

Power Convortor Modnl ing Ohjoctiv os 

1. Accuracy of terminal characteristics 

2. Accuracy of AC model for stability and dynamic performance 

3. Simplified internal working model 

4. Capable of DC, AC, and transient analysis 

5. Time domain model for large signal simulation 

6. Frequency domain model for small signal analysis 

7. Capable of predicting system interactions such as input filter 
Interactions 

8. Equivalent circuit model and transfer function model 

9. Model adaptability to canned circuit analysis programs such 
as SPICE 2, SCEPTRE, ICAP, etc. 

The most applicable modeling and analysis technique depends on the 
following considerations: 

1. The analysis objectives: worst-case analysis; DC analysis; AC 

analysis; transient analysis 

2. Accuracy required 

3. Type of control circuit used 

4. Nature of disturbance 

5. User's background 

6. Capability of the host computer 

P ropo sed swi tching Co n v erte r Model in^ Schemes 

1. DC and small signal AC model. 

Power stage: discrete-average model 

Duty cycle modu1.Htor: describing function 

I t'niv |iiu>essor; • ^ air ler t ;rU.t uiti ''lodel 

2. large '.ignal model . 

P'lwer stage: a) discrete-average model for cost-effective 

analysis and simulation 

b) diM.rete time d'eiiiuin simulation 

i) ditrU lirciiit si'i ilal ion programs 
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Duty cycle modulator: direct simulatioii uainn equivalent circuit 

model or nonlinear differential equivalent 

Analog error proces- direi:l, simulation using equivalent circuit 
sor: model or differential equations 

The particular chosen modnling tecliniqiies should be simple and effective 
for largo system simulation. On the other hand, it should bn accurate enough 
for detail study and trouble-shooting. Tliose apparently conflicting objec- 
tives require inultiplo converter models to be ostablirdied and a particular 
one is selected for a specific application, 

Dev el 0 pine n t (• f f o r t s 

The development efforts needed to establish proper modeling tools for 
the proposed modeling schemes are summarized as follows: 

1. DC and small signal model 

Power stage model : 

t The discrete-average technique is currently developed under the 
support of NASA Lewis Research Center, Grant No. NSG-3724. 

t Discrete-average model for three basic converters; buck, boost, 
and buck/boost, have been demonstrated. 

• The discrete-average model should be represented in the follow- 
ing forms for generality: 

- state space model 

- transfer function model 

- equivalent circuit model 

• The Miodel ing scheme should be extended to other convener type*"- 
to estatihsh a componetit (lonverter) library 

Duty I yi ie itiodul a t nr model . 

• Ihe accu.acy ul Uie duly Lyck iiiodulalnr model is 1 in ed to 
only low modulation frequencies using the describing function 
technique. 

• The inaccuracy of the pulse modulator model is particularly 
pronounced for high performance converters employing multi- 
loop control . 

• An improved modeling sclu'iiie is desirable. 

2. Large signal model 

Power stage model : 

• Discrete-average model - advantages: 

A simplifind inti'rnal working model; the imv|el is simple, easy 
to u'.(> , and 1 11'. I -(■ 1 1 1'l t. i VO . 


— r. 
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• D1scrot(.*"av(M’dO() iiiodol - dlsndvtintiioo'i : 

i III! iiiodol is midbli! 1.0 liuridlo abiUH’iiuil toiivor l.cr opoidiions 
Kui;h AS poAk curroni: protection mode of opnratlon; f.Aluratlon 
of iiiAiinotIc compononts; saturation of op-amp In Uio foodhack 
control loops. 


• Dlscroto tliiio domain analysis A simulation - advantaflos: 

fxact duplication of circuit holiavlor capablu of Incorporating 
all system noni Ini'arltlos. 


Accurate and cost-offuctivo. 


Discrete time domain analysis & simulation - disadvantages; 

Relatively Ineffective for large scale system simulation. 

Direct circuit analysis & simulation - advantages: 

Easy to Implement, flexible and versatile. SPICE 2 and SCEPTRE 
are available CAD programs for converter simulations. 

Direct circuit analysis & simulation - disadvantages: 

Time consuming for Implementation and numerical integrations. 


? . 1 . 5 C a bj i ng and Distrib u 1 1 o n 

Table 2-9 summarizes the capabilities of the distributlon/cabi ing 
modeling techniques with respect to various attributes. Table 2-10 
summarizes the existing model uncertainties/inadequacies. Cabling and 
distribution optimization models fall Into two classes: 

1. Mass optimization 

2. Cost optimization 

In the first case the optimum distribution weight is determined as a 
function of conductor mass and conductivity, the specific mass of the 
source and storage elements. The existing models reviewed fall to 
account for fuse, switch and connector losses, which in some cases could 
override conductor losses. 

In the second case the optimum conductor cross sectional area is 
determined as a function of conductor, energy storage, solar array, 
Ihrniicl control, and power conversion costs. In general, the optimum 
lost abling is heavier Ihiin the iiiin iinuiii mass node!. The cost model 
has no provision for insuring a cable at least as large as the mi-. .mum 
Ilia S'.) .node I . 
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DISTRIBUTION 

MODEL 

ATTRIBUTES 


• MODELS EOR I, 8 

• MODELS FOR 2,3,8 

• MODELS FOR 4,5 

• MODELS FOR 6,7 

• MODEL FOR 9 
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Component Kodpi Capabu ity Summary 
ni'iTPinmioN/rAiH inp 



• GENERALLY ENTERED AS A WORST CASE RESISTANCE. TESTING WY BE REQUIRED TO 

OPERATIiNG CITARACTERISTIC8 OVER A WIDE RAfftSE OF ELECTRICAL AND 
THERMAL CONDITIONS 

• LINE INDUCTANCE HAS BEEN MODELED, NO LOOP INDUQANCE MODELS WERE FOUND, 

“ PRECISE INPUT DATA IS REQUIRED FOR TRANSIENT ANALYSIS 

• THESE HAVE BEEN MODELED ON A SMALL SCALE 

- PRECISE INPUT DATA IS REQUIRED FOR TRANSIENT ANALYSIS 

• NO COMPLETE MODEL EX' STS, SOME MODELS SIZE BASED ON MINIMUM MASS, 

OTftRS SIZED ON MINIMUM COST 


Table 2-9 
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Distribution Mooet. 
Uncertainty/Inade-.qijacy 




DIFFICULTIES OF ESTIMATING DISTRIBUTED 

MuruAL coupling/cable routing effects 


INDUaANCES; CAPACITANCES DUE 


TO 


• VARIABILITY OF CONNEaOR, SWITCH CONTACT RESISTANCES .WITH AGE 

• SENSITIVITY OF HYBRIDS AND l/c's TO TRANSIENTS 

MAKES TRANSIENT PREDICTION AND ANALYSIS DESIRABLE 

- MAKES DIRECT TESTING UNDESIRABLE DUE TTD DANGER OF 
PART OVERSTRESS 


j 


Table 2-10 
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2,1,6 Mpdpl 1 ng Of Cpinpl of o Powpr 

Sovoral qonorrtl cirtc.r.os of powor finbsy'ifom (iiodol Inq proprdiiui sppoar* 
In tlio 1 If.oraturo, iho'.o aro; 

1. (lorioralUfKl rlnutlf oikI uyntoiiifi analy.is inodolr., 

)•' 2. Dodicatod p'leudn'irtiphlcal iiiodols of 'ipoflflf. ‘iulMiysfoiTtfi. 

3, Ilyhold »’oii)p(jtaf.1on modoln. 

:» Tablo ?"li suiiinio*'’ Uos tlio capabi 1 of ‘jpoolflc models sl.udlod, 

h> and Tublo 2-12 simiiiuirU’os tho adopuucy/1 Imitations of power subsystem 

model s . 

2.1.6,! Genera ] j ^ ed Cpmjju ta y p iia 1_ Mo^^ 1 s 

Ihe generalised computational models, as a class, are represented by 
programs such as SPICE, SCEPTRE, ECAP, ICAP, TESS, and a number of others. 
Variations and extensions to these programs are also reported. 

All have similarities In that Lh(?y offer a set of typical component 
models which are assembled Into a network of interconnecting nodes. Each 
individual component model must be reducible to one or more first-order 
differential equations which are solved by matrix arithmetic to find a 
solution at each point In time. 

Many of these programs have the ability to accept parametric data In 
the form of equations or functions, piecewise linear parametric data, or 
as user-supplied FORTRAN routines. The supplied equations and routines 
must meet the criterion of being reducible to differential equation form. 

When used with accurate models of individual components, these oro- 
grams are generally capable of performing non-linear DC and transient 
analyses, and AC analyses if the AC signal Is small enough to permit tho 
assumption of linearity within the range of interest. This class of 
general k'ed programs shares certain limitations: 

1. Witli increasing numbers of components, the memory requirements 
Increases dramatically. By use of dynamic memory management It 
Is possible to permit an unlimited problem size, but only at 
the cost of an Increase in computation time. 

2, Aside from the memory management problem, an increase In problem 
size results in a geometric increase in computation time. This 
makes the analysis of large, complex systems expensive. 


oniG!N;a, FACir? n 
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POWER SUBSYSTEMS 
REFERENCE 





'W) 






M/ 


14 Mdx 1 ml 2 lng S/C 
Energy 
Ostwald 


l 7 t>'Plorer Energy 
Balance Program 
Broderick 


22 ?wr. Systems, GEO 
Blllerbeck 


26tnergy Balance 
Program 
NA^.A/GSrC 


Note: References are keyed 

to AppenfHx A 


A 



/ ^ I 

'iS h 


SELECTED POWFR SIIRSYSTEM MnDEl S CAPABII ITY SUMMARY 


0-Has capability 
x-LlmIted capahil tty 
Rlank-Nn rap.thilitv 


TABLE 2-n 










NASA CR^16682Q 
TRW .18651-0001 


Adequacy/Limi TAT ioNs op 
Power Subsystem Models 


• ability to do small signal AC ANALYSIS 

• NO ABILITY TO DO URGE SIGNAL AC ANALYSIS 

f LIMITED ABILITY OF MODEL TO FUNCTION IN TRANSITION 
ORBITS OR ELLIPTICAL ORBITS 

• ONLY ONE MODEL WAS FOUND TO HAVE RESTART CAPABILITY 

• SOLAR ARRAY MODELS UCK PUSMA INTERACTION DATA 

f SOME MODELS CAN BE USED FOR SIZING^ HOWEVER IT IS 
GENERALLY COST EFFECTIVE TO USE A SIMPLER MODEL 

e MODEL VALIDATION DATA IS UCKING . THIS IS REQUIRED 
IF THE MODEL IS TO BE USED AS A SYSTEM VERIFICATION 
TOOL 


Table 2-12 
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By permitting the entry of data In the form of functions and subroutines, 
the more modern versions permit the summarkatlon of the performance of large 
numbers of components by a sitigle (or a small set of) transfer functions, thus 
reducing the problem size for large system to more manageable proportions. 
However, the program then becomes more difficult to run, for a suitable transfer 
function for each "box" or circuit must be found, which adequately describes 
Its operation under all conditions within the range of Interest. For some 
components this can become a formidable task. 

2. 1.6. 2 Dedicated Pseudographical Model g 

Dedicated models of Individual power systems have been used successfully 
In many cases. In these the logic and properties are programmed into the 
computer model and may be In a variety of forms, including logical equations, 
subroutines, and piecewise linear approximations. One Is not limited to any 
specific form, nor Is the ability to form a first-order derivative essential, 
although differentiabll ity may be a convenience in the use of some numerical 
solution routines. Because of the lack of adequate mathematical descriptions 
of some power subsystem components such as batteries, the solution method 
used will usually take a pseudographical form, i.e., the computer will find 
the solution by finding the intersection of two tabular functions. 

This type of model is capable only of DC analyses, although if heat 
generation is considered, and if the thermal time constants are much longer 
than the electrical time constants (which they Invariably are), the system 
will follow thermal transients in an electrical system very well. 

Program accuracy is limited by the accuracy of the model logic, the 

input data, and the accuracy of the Interpolation or spline routines used 
for finding the solution. 

Execution speed and cost of such programs tends to be much less than 
that for the generalized circuit and system routines for large systems, and 
the accuracy greater, primarily because the design flexibility of the indi- 
vidual component models Is greater. 
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2 . 1 . 6 . 3 H>' b r1 d Mp d o l r> 

In hybrid modeling the mathematical transfer functions used to simulate 
the porfornianco of circuits or boxes are replaced by analog computer models 
having the same function. The analog-modeled transfer function outputs are 
used In the same way as the mathematical transfer functions In a generalized 
type of program. The result in theory is to speed up the computation (as 
compared with all-up part models). The method has an additional degree of 
flexibility. By placing the parameters of each analog model under the 
control of the digital program, it Is possible to vary the transfer functions 
of several such models simultaneously so as to converge toward a desired 
solution. At this stage, it is not clear how far this method has been 
developed. 


indu stry S urvey 

In order to determine the actual usage of models, an industry survey 
was conducted. Because of the nature of the subject and the limited time, 
not all companies were contacted, and some of the data is necessarily brief. 
However, the message is clear that no overall comprehensive model is in 
existence, and no effort is underway to develop such a model. 

Aerospace Corpor at i o n 
• General purpose solar array model 


- Input: 

- Output: 

- Limitations: 


Array geometry, sun angles, radiation dose rate, solar 
cell parameters, temperature characteristics, mission 
1 ength . 

I-V table, and solar array I-V curve plotted, 

DC model. Array temperature and distance from sun 
must be determined by user as an input. 


Power system simulation program 

Variation of TRW's PSIM program. No transient analysis capability. 
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BOEING Am’p t; n 

• Sliiipl 1 find s1/1tu| program 

- Battery and solar array characteristics have to be changed to reflect 
degradation and aging. 

t Performance verification program 

" Input: Orbit parameters, load profile, attitude pointing pro- 

file, solar array pointing profile, max battery DoO, 
line losses. 

Output: Thermal data. Including transients, DC analysis, voltage, 

current, state-of-charge, system sizing requirements. 

- no electrical transient capability 

- developed for direct energy transfer and centralized 
regulator systems 

Ford Aerospace 

• Improved SPICE (1 SPICE), ECAP 

- Developing detail models for standard converter design. 

Lockheed Missiles and S pace 

• Sizing model 

- Proposals, initial concepts 

• Performance Verification 

- On-going spacecraft 

- Mission specific 

- No transient capability 

- Verified by test 
Gen eral Elect ric 

• Energy balance program - low earth orbit 

- Steady-state DC 

- No transient analysis 
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Marti n Mar iett a ( Qonvor) 

• Photovoltaic systoni tost prototype model 

- General purpose terrestrial or space, with or without battery 

- No transient analysis 

- Verified by tests 

• Developing new version of SPICE 

• 5KYLAB computer program 
RCA 

• Energy balance computer program 

- Multiple orbit capability 

- Computer system 

- Linearized electronic component models 

- No transient capability 

- No interactive thermal model 
TRW 

• SOLAR 

- Predicts solar array output I-V curves 

- Accounts for degradation factors 

- Accurate, easy to use 

- Interactive, rapid turn around 

- Validated by flight data 

- Requires multiple runs for complex array geometry 

- Requires service program to input tape for power subsystem program 
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TRW (Cont'd) 

• SOL GRAD 

- Same as SOLAR except 

-- two-dimensional thermal gradient model 

— provides tape fur power subsystem program 
-> more complex data input required 

• HOTSPT 

- Predicts heat generated 1n solar cells as a function of shadowing or 
breakage, and the resultant solar array 1-V curves 

— not Integrated into power subsystem model 

— assumes all cells have same reverse characteristics 

• BATMODL 

- Predicts current, voltage, heat generation, Internal pressure, 
soc of battery cells in low earth orbit 

- Not fully tested or validated 

• DISTRIB 

- Determines minimum mass cable for a given set of loads and power 
sniirre mass per wa'f't. 

• PSIM 

- Energy balance/ power system performance computer program 

- Uses solar array tape from SOLGRAD 

- Models complete power system with load and thermal interfaces 

- Restart capability 

- Voltage, current, heat generation, battery SOC 

- Secondary distribution networks not modeled 

- No AC or transient capability 

- Used for the following spacecraft: HEAO, FSC, DSCS II, DSP 
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2,1,8 Ideal Power. Systqnis Model Roi|u1ri!iiientu 

The first requiromnnt of an Ideal model Is that It must he user-orlonted. 
The model must he capable of performlnq conflquratlon trade studies to de- 
termine tlie optimum power system for a particular spacecraft. It must bo 
capable of analyzinq the available options and quantifying the results. 
Optimization may be on either cost or mass, or on some other parameter, such 
as minimum solar array size, or minimum radiator size. This flexibility has 
to be comprehensive enough to minimize changes downstream in a project wheeo 
It can become very costly. An output of the trade studies Is the Initial 
power system sizing. 

- solar array power and area 

- energy storage mass and dissipation 

- distribution system mass 

- regulator mass and dissipation 

The model must also be capable of identifying design and performance inade- 
quacies as the subsystem requirements evolve. Early prediction of potential 
problems will minimize design changes. 

To insure development of stable and reliable power systems the model 
must predict system stability margins and be accurate enough to allow sensi- 
tivity analyses to be performed. 

The above requirements can only be met by a very sophisticated computer 
program with very accurate models of all of the power system components and 
loads. This accuracy is not required for initial trade studies and Initial 
power system sizing. A much simpler program could do the job, and it appears 
that It would be used for specific applications. An example would be to have 
separate programs to compute: 

1. Initial trades and sizing. 

2. Nominal performance predictions. 

3. Off nominal or design limit verification. Failure analysis, 
and performance in a degraded mode. 

4. Transient analyses - system stability 

These could all be contained on a tape and the desired function and models 
called by command. 
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3.0 COMPRI.fIKNSIVI POWIR SlIllSYSTLM MOOU-ING APPROACH 

Roc]iii ronton Is 

Tlio usos, and thoroforo tho roqulroniont s , of n pnwor ‘.iihsystnm model 
vary with tho phase of the spaoeoraft devolopnienl prooram in which they are 
applied, fiijure 3-1 outlines the work flow plan followed for this task. 

Protjrani Phase A: Concept definition and selection 

In this phase of the program neither the requirements of the power 
subsystem, nor its configuration, are well-defined. The effort of the con- 
tractor is to define the requirements and characteristics of a power subsystem 
and to select a configuration which will meet those requirements. The 
customer has the task of evaluating the selected configuration to determine 
whether or not it is best suited to the requirements. In this phase, only a 
sizing or synthesis model can be used effectively. Performance models cannot 
be used effectively until the design has been advanced to the point at which 
a reasonably accurate model can be constructed. 

^se A Requ irements 

The Phase A model must perform at least the following tasks: ' 

• Provide an estimate of the relative cost, weight, volume, solar 
array area, and other optimization criteria for a large set of 
possible power subsystem configurations. 

• Provide for each power subsystem candidate configurations, inter- 
face information necessary for the definition of the requirements 
placed upon other spacecraft subsystems by the power subsystem. 

(For example, heat dissipation, mass distribution, torque induced 
by rotating machinery, etc.) 

Many such models already exist. Most are proprietary, and documentation tends 
to be relatively poor. This type of model is a relatively simple one to con- 
struct. and represents little technical challenge. Its usefulness lies 
primarily in its ability to 

• Provide rapid response 

• Minimize Phase A slady costs 

• Avoid overlooking some important factors 

• Provide a vehicle for lustoiiu'r evaluation of proposed 
con t i«|urat ions . 
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rrPirfll'I -Piwri/* O. Pro.1.1n|1,nary , Dor.1(jn 

In Phafio n, thfl dosiqn concopl liaf. boon i,(!l(H'locl fiiul tho niajiir lonl.rat.to* 
nffort Is dlroctod toward tho dovolopinont of Imploiiiontat'lon coiicopts. bnno- 
ral circuit doslcm typos aro roiHihod out. hut coiiiponont valuns have not all 
boon soloctod, and tho dofinitlon of tho powor rotiulroinonts. wh11(> Iniprovod. 

Is still not firm. Uttlo or >\oth1n<j Is usually known about tho individual 
loads other than thoir approxliiiato power consumption and tho required Input 
voltages. Circuit definition and power distribution network definition has 
not yet reached the point where either AC or transient modeling Is possible 
In a meaningful way. Only If the power .subsystem consists of existing hard- 
ware which has been extensively characterised can the AC and transient modeling 
be done at this stage. Even then, the spacecraft-unique distribution system 
win be unavailable for use In an all-up model. A DC model can be extremely 
useful In Phase B for tho definition of many of the power subsystem character- 
istics. 


Mo del Requirement s 

The Phase B DC model is used primarily to verify the power and energy balance 
characteristics of the power subsystem, and to provide assurance that there 
is sufficient power and stored energy to meet mission requirements, given 
the assumed characteristics of the power subsystem components. The require- 
ments of such a model aro; 

• Predict voltage, current and heat output of each power subsystem 
component (including batteries). 

• Predict the voltage, current input required for each component 
to deliver the required output. 

• Predict voltage drops in the interconnecting harness and power 
distribution harness. 

• Predict energy content of energy storage devices, such as batteries 

• Simulate the electrical and thermal response to 

- load changes 

- loss and return of solar array Illumination 


comma nds 



NASA CR-1 66820 
TRW 38661^0001 


- I'liv ironniontid chiuiqos (‘iuch ns ‘‘lOlnr Input, dlhoclo, 
r,/c orlnntntlon) 

t.pai r. cimriourat Ion dianrp's (lintoro aiul aftor 
(Joployiiinnt) 

" pnxTIct.thlo (’(rllur<;s. 

Prod'lctlonr. an> niddo on a functlnn oT tlino. Tho t.otiil tliiic npan of 
itil'.orofit 1'j one of multiple orhitn. 

d.. F- Onto 1 1 od Den 1 yn 

During Phase C the designs arc? completed, All parts have been defined, and 
both electrical circuit dcjsigns and physical layouts have been complotod. 
Engineering models have been constructed, and some measurements have been 
made for the purpose of characterizing their performance. These measure- 
ments provide a more accurate estimate of power subsystem component DC 
performance, permitting Further DC performance predictions for the entire 
power subsystem at considerably Improved accuracy. 

Also during Phase C, sufficient detail is generated to permit the construc- 
tion of AC and transient models of the components of the power subsystem. 

Early in this phase, such models arc useful In performing trade studies on 
alternative detailed circuit designs. Toward the end of this phase, AC and 
transient models are used to analyze the interaction between components of 
the power subsystem. 

P hase C Mode l Requ irem ents 

Requirements are the same as for Phase B, except that the accuracy 
with which some of the power subsystem components have been characterized is 
improved as a result of measurements. 

.F^Silnaj_ AC The primary use of small-signal AC models is a de- 

termination of component and/or subsystem stability. Suc:h a program should 
be capable of performing the following predictions: 

• Determination of component stability. Such predictions are usually 
in the frequency domain in tcriiis ot g,iin and phase margins. As 
such, they are independent of time. 

• Individual control loop characteristics, stability, gain, phase 
margins . 
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• Input-output. Iinpndnnr.o nr, n funct ion of froquoncy. 

• Smqll-slqndl trnnsinnt rorponso, 

• Cflpflblllty of Incluclinq lonipoiumt intoroct loni. , ptjrticuinriy 
thor .0 of input mid output flltors, 

j-fft’SPrSJijnal fron'. lent Mod(?1 1 ; ! L lko thn DC niodol, the 1nrqn-r,1(|rwi1 tranr>1ont 
iiiodol 1r, non°l Inonc, and oporato!. In tho t1nio doiiiirln. Howovor, for ocnnoiirlr. 
n'cisoivi, tho trails I out iiiodnl Is lliiritod In Its tliiio span l.o porlods of vary 

short duration; of tin; orclor of mill Isutomls. Tlio uso of such translimt 

iiioduls Is to predict the voltarif) and currcuU translGnts which occur as a 
result of an outside stimulus or Internal failure of a component. Tho result- 
ant prediction Is then axaiirinod to determine whether or not It will have a 
deleterious effect upon the operation of the component or other components of 
the system. The requirements of such large-signal transient models are as 
follows : 

• Predict the voltage and current transients at the Input 
and/or output of a power .subsystem component (or set of 
components) ; 

as a result of an external stimulus, 

as a result of an internal failure. 

Program Phase D - Hardware Construction and Spa re craft Integration 

In this phase, flight hardware Is constructed and the complete spacecraft is 
assembled. Measurements are made characterizing the performance of the 
Individual components of the system. In smaller spacecraft, an integrated 
subsystem test is usually performed to determine the overall performance of 
the system. However, with increasing spacecraft size and Increasing power 
levels, the difficulties inherent 1n such integrated subsystem tests multiply. 
No new kind of subsystem model is needed in Phase D. The same DC, small- 
signal AC, and transient models are used. However, more accurate results 
should be attainable because the measured characteristics of individual flight 
hardware components are used as Inputs in place of the earlier estimates of 
component performance. 
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Program Phano E - Flight Qpnrati onn 

In thin phano tho npacncraft In placod into uporaMonal orbit (or flight 
path) and In opnratod, Tho iinon of powor nubnyntoin modoln In thin phano 
of npacocraft oporatlon aroj 

• Falluro dlofirtontlcn (ind prodlrt1onn of nff-noinlnal oporatlon 

• In tho onno of In-flight inaliitalnaldo npacooraft , prodictlonn 
of ohangon Is tho porforiiiaiion of tho powor niilmyntoni an a 
ror.ult of rhangon 1ii l;lio npanooraft powi.'r domand or powor •iOb“ 
nyr.toni conligiirntlon, 

• S.ynfoiit hortlth niatii‘i/<,onipotiont por lorinanoo. 

Although no now typon of programn arc; roc|u1rnd to accoinpl 1t.h thono a linn, 
additional roc|Ulreinonts may bo placod upon oach of tho program typos to 
Include the following capabll 1t1o;; : 

• Model a fairly extensive set of plausible failure inodes. 

• Follow changes in spacecraft load and power subsystem configura- 
tion. 

Figure 3-2 shows a block diagram for a power system model. 

^ ^ J A ttr ibute Re qui remen ts 

In the preceding paragraphs, we have defined the things the models 
must do. In the following are defined certain essential attributes or 
characteristics which all of the models must have in common to create a 
comprehensive modeling and analysis capability: 

t Commonality and compatibility: Each of the types of models must 

maintain a commonality of reference wHIi all of tho others; I.e., 
they must all predict the same attributes and performance for the 
same subsystem. 

• Modularity: Each of the components of the power subsystem should 

be modeled as an independent module having the following attributes: 

Capable of being operated independently as a conipcnent 
model 

Capable of being integrated into a powor subsystem 
model 

Capable of replacing or being replaced by an alternative 
model having different input requirements 

Data base - each module must be provided with an independent 
data ba'..t> specifically suited to its own needs. 
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• Efficient use of both core memory and computation time. 

• Vori f iabi 1 i ty - models should be verifiable as Independent modules 
and as a complete integrated power subsystem. Verification data 
base must be provided independent of the model input data base. 

• Operational Simplicity - The models must be designed for use by a 
working power subsystem engineer whose understanding of the program 
and computational facility are limited. Clearly stated, complete 
user documentation is essential to the successful use of such a 
complex set of models. 

3 . 2 Single Comprehensive Program vs Independent Programs 

The need for four types of programs has been identified, each of 
which has different uses, and uses different methods to accomplish its 
objectives. It is possible, however, to design a single program containing 
the four programs which will accomplish all of the required objectives 
without severe penalties in either computation time or memory requirements 
by the use of overlays and segmentation methods. By this means, only those 
parts of the program which are being used are resident in memory at execu- 
tion time. 

3.3 Proposed Modeling Concept 

The comprehensive Power Subsystem Model is envisioned as a set of 
four programs, either entirely separate, or combined into a larger segmented 
program, the desired part of which is called into core memory at execution 
time. Each program is equipped with a set of component models. These are 
selected by the user and are interconnected by the program to form an inte- 
grated power subsystem model . Figure 3-3 shows an overview of program 
structure. The operation of the program is as follows: 

• Based upon the case instructions and input data, the executive 
selects the desired program segments and the specific power 
subsystem component models required for the run, and copies them 
from the disc library. 

• An object program is generated, containing only those program 
segments and component models necessary for the specific run, 
and loaded into memory for execution. 
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• lilt’ il.ila niiUl.iycM', iiiuhT ('xoi iit Ivi' ilirtH i ioii, ‘-(’Ici t s l.hc i'(n|'i Irod 
in(Hit il.il.i lor iMch oT I ho iioinpotuMil iiiodols, ro.uls I ho daf.i from 
disc packs into a loiiiiiioii slorayo area, and (orms an index of daia 
local ion. flopondifui upon Mio si.’o of dal.a hanks, randoiii-accoss 
dii.c opiM’ations may h(’ sidisl i1 nl cd for in-cort' data storayo, 

• Ouriny proyrum oxocmion, c.ompononl model oiilput. data are yono- 
ratod, stored in the computer memory, and indexed under the 
lontrol of tile data manayer. Tlu\y may subsequently be written 
to disc or tape. 

• At appropriate points in proyram execution, tlie output data are 
retrieved from storage and printed or plotted. 

A program structure of this kind represents a compromise between minimizing 
computer RAM memory size and program execution costs. As the data bank 
sizing and data manager execution costs become clearer, an optimum memory 
field will b ecome d e f i na b1 e . 

d . 3 , 1 Computer 1. a np uaje 

Prior to coding, consideration will be given to the selection and use 
of computer language. However, unless significant advantages are found in 
one of the other scientific languages, l-ORTRAN 77 will be used. It is felt 
tliat this language is the richest in available internal functions, and the 
most commonly understood of all languages in use in the scientific and 
eng i neei' i ng commun i ty . 

3 . 3 . c 0 i sp 1 ay forma t s 

Considerat ion vfill also be given to llte development of a program control 
format based upon MINU selection of program functions. This is desirable on 
systems containing CRT displays, and the potential advantages of such a system 
in terms of user simplitity and compat inil ity are great enough to deserve 
considerat ion. 

3.3. ,t Program Di've I opmen I 

figure .1-4 shows a sugiiesled sclu'dule for model and ilala base' develop- 
ment. This schedule is flexible, and is c.ipalih' ol rc'sponding to variatii>ns in 
funding availability, NASA priority of inlero'.! ami otiier f.ictors. It stiould 
however follow tin' gi'iieritl flow ol I inure .i-'.*, and is sub.iect to certain 
constraints upon the (trd.'r in whicli the variou'. pli.e.e'. ol I he efiort are 
carried out. 
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Driver Progr am Dovo lo pnient 

It. will be necessary to complete a detailed definition of the driver program, 
data base structure, and data manager before beginning coding of the indi- 
vidual models of power subsystem components. Failure to do this could result 
in severe incompatibility between component models and the driver which 
operates them. Coding of the driver can be delayed, to some extent, but 
should have been completed and debugged prior to their first attempt at 
running an all-up power subsystem model test. 

Model Developm ent 

Models may be divided into two categories: existing models which require 

only adaptation to the program structure, and new or improved models which 
require development. New model development follows the flow of Figures 3-6 
through 3-10, proceeding through trade studies to select the approach to be 
used, development of the mathematical or logical algorithm, and coding. 

Trade studies and algorithm development can be done prior to or concurrently 
with the program structure development. Coding must await the development 
of program and data bus structure before it can be accomplished efficiently. 
Component models are independent and modular, so that updating a component 
model does not require changing the code of the executive program. 

Data Base Development 

For those models whose concepts are well understood, a data acquisition 
program can be initiated independently of the development of the final model 
code. However, where the model algorithms have yet to be developed, the 
form of the data required, and in some cases the kind of data required, to 
provide an input to the model are unknown or ill-defined. For these cases, 
it will be necessary to defer the development of a data acquisition program 
until the algorithm to be used has been defined. Data base acquisition 
requirements are shown in Tables 3-1 through 3-5. 
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A Suggested Battery ^^0DEL Development Approach 


NASA CR-166820 
TRW 38651-0001 


ORic.irr‘\ [ •; 

OF POCIJ QUALUY 



FIGURE 3-6 






















A Suggested Solar Array Model Development Approach 
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FIGURE 3-8 














A Suggested Shunt Regulator Model 
Development Approach 
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FIGURE 3-10 
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Data Base Acquisition - Batteries 


• SUGGESTED ORDER OF TESTING 

• EFFICIENCY/SOC 

• HEAT generation/pressure 

• CYCLING PULSED LOAD/CHARGING 
- AC MODELING 

• AGE-LIFE 

• HIGH VERSUS LOW VOLTAGE 

• CELL SIZE 

• high/low TEMPERATURE 

• VENDOR-VENDOR VARIATIONS 

• SUGGESTED TESTING OF BATTERY TYPES 

• NiCd 

• N1H2 

• AgHg 

• OTHER 


TABU' 3-1 
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Data Base Acquisition Power Distribution 


• SUGGESTED ORDER OF TESTING 

• NEW DEVICES 

— SOLID STATE SWITCHES 
— SEMICONDUCTORS 
— HYBRIDS 

— INTEGRATED CIRCUITS 
— ON-GOING TECHNOLOGY 


TABLE 3-2 
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Data Base Acquisition - Solar Arrays 


SUGGESTED ORDER OF TESTING 

• RADIATION DAMAGE^ PROTON DOMINATED ORBITS 

t COVER GLASS AGING/ ADHESIVE AGING 

t TRANSIENT RESPONSE 

f THERMAL CYCLING BEYOND 5 YEARS (> 30/000 CYCLES) 

• SPACE PLASMA INTERACTION 

CONCENTRATOR ARRAYS 

• SPECTRAL REFLECTANCE OF SPECIFIC COATINGS 

• MIS-ALIGNMENT AFFECTS 

' EX THERMAL EFFECTS ON CONES/ CELLS AND STRUCTURE 

al-optical properties of surfaces 


TABLE 3-3 
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Data Base Acquisition -• Shunt Regulators 


9 SUGGESTED TESTING 

9 NEW DEVICES - ELEMENT PARAMETER LIBRARY 
— SOLID STATE SWITCHES 
-- MOSFETS 
— HYBRIDS 
— CONTROL CIRCUITS 
— ON-GOING TECHNOLOGY 

• DEVELOP ELEMEI^T PARAMETER LIBRARY 

• FOR DIGITAL SHUNTS^ DEVELOP METHODS IN TIME DOMAIN 
USING SAMPLE DATA TECHNIQUES 

TABLE 3-4 
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Data Rase Acquisition 
DC-DC Converters 


• SUGGESTED ORDER OF TESTING 

• LARGE SIGNAL AC, TRANSIENT RESPONSE 

• TRANSFORMER 

- COUPLING 

- LEAKAGE 

- SATURATION 

• DEVELOP ELEMENT PARAMETER LIBRARY 

• ON-GOING TECHNOLOGY 


TABU 3-S 
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i».u coNaiisioNs 

Ki(> I It oratiiro ‘.('.u’i h/liulusi ry survoy l om inntul th.it no t om(>ri'lu'iis ivo 
powoi* ‘iiilisysloin inoilol ox I si ;.. 

It In (lilt 1i i[M toil thill .1 sol of rotio rmul.iinoivliil t.v|.io;i ol iiuiih'l;. .iro 
i'iH|Uli'ocl, o.u'h of which poi’I'oniis .i dlffori'iit , onnoiit i.il I, ink. liuioUior, 
this pi’iniiMiii st'l. lUM'iiilts coiiipi'i'hoiin ivo iiioilolliu) of a powoi' siihsystoiii. 

1. A powoi* siihsyntoiii ;il/iii«| ami nynthoni;; proiji'aiii, oap.iblo of 
ostliii.it I lip cost, mass, voliiiiio, ai'o.i, and othor attrlhiitos 
of a sliipio-point doslpii. This would ho usod only dui’liuj 
tho concoptuul dositjn phaso of a spacoci'aft poopi'am. 

1*. A lie iiiodol of tho powoi* siihsystoiii and Us intorfaclnp suh- 
systoiiis. This is usod duoinp phasos B, C, D, and f ol tho 
spacocraft propoaiii for studios of powoi* consumption, ooonjy 
halanco, hoat ponooiiHon and absonption, rosponsos of tho 
suhsystoni to oiivironmontal chanpos, and prodictlon of 
si oady-stato voltapos and cui*i*ont thoouphout tho suhsystoni. 

3. A small-sipnal AC model of tho powoi* suhsystoni ami its 
compononts. This Is usod in phasos C, l\ and I for tho 
pui'posos of ostimatinp suhsystoni and intorcomponont stability, 
bus impodanco, and for dotonnlnation of rosponsos to small- 
sipnal trails I outs. 

4. A larpo-sipnal transiont modol usod during phasos C, 0, and 

I- of tho program for tho purpose of dotoriiilninp tho response 
of tho suhsystoni to largo transients such as state changos 
and faults. 
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DC-DC Converters 

A-10 


The following abbreviations are used in the references: 

lECEC, Intersociety Energy Conversion Engineering Conference 
IEEE, Institute of Electrical and Electronic Engineers 
PCSC, Power Conditioning Specialist Conference 
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Da ta Base : 

Battery models used IRAD and 
Industry data. Same for solar 
array/ 




TECHNOLOGY EVA - TION SUMMARY - TELECOf 



Constraints were time and costs 
therefore, simple model. 

Limitations : 

No signal ac 

No transient analysis 


Outputs : 

To be used for 6 major subroutines: 
Attitude control, solar array. 
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